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HERMITE-PADE APPROXIMANTS FOR A PAIR OF CAUCHY TRANSFORMS 
WITH OVERLAPPING SYMMETRIC SUPPORTS 


ALEXANDER I. APTEKAREV, WALTER VAN ASSCHE, AND MAXIM L. YATTSELEV 

Abstract. Hermite-Pade approximants of type II are vectors of rational functions with com¬ 
mon denominator that interpolate a given vector of power series at infinity with maximal 
order. We are interested in the situation when the approximated vector is given by a pair 
of Cauchy transforms of smooth complex measures supported on the real line. The conver¬ 
gence properties of the approximants are rather well understood when the supports consist 
of two disjoint intervals (Angelesco systems) or two intervals that coincide under the condi¬ 
tion that the ratio of the measures is a restriction of the Cauchy transform of a third measure 
(Nikishin systems). In this work we consider the case where the supports form two overlap¬ 
ping intervals (in a symmetric way) and the ratio of the measures extends to a holomorphic 
function in a region that depends on the size of the overlap. We derive Szego-type formulae 
for the asymptotics of the approximants, identify the convergence and divergence domains 
(the divergence domains appear for Angelesco systems but are not present for Nikishin sys¬ 
tems), and show the presence of overinterpolation (a feature peculiar for Nikishin systems 
but not for Angelesco systems). Our analysis is based on a Riemann-Hilbert problem for 
multiple orthogonal polynomials (the common denominator). 
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1 INTRODUCTION 


For p £ N, let f = (t-\,, f p ) be a vector of holomorphic germs at infinity. Given a multi¬ 
index ft = (ni,..., Up) £ N p , the type II Hermite-Pade approximant to f corresponding to 
n, 




(1.1) 


7T* = 


n [ } ] ■= UL. 
" ' Qrt' 


K”.O- 

is a vector of rational functions with common denominator Q b satisfying 
deg(Qrt) < |ft| =m -1-+ Tip 

(1.2) 


(z) := (Qnfi-Pj ) )(z)= 0 (z-^- 1 ) 


R^iz] := 


as 


00 


for each i £ { 1 ,... ,p}. We shall not deal with type I Hermite-Pade approximants in this 
work and therefore henceforth we will drop the “type II" modifier. Such approximants 
were introduced by Hermite [23] for the vector of exponentials ( 1 ,e z ,..., ef p_1 iz ), with 
the interpolation taking place at the origin rather than at infinity, as a tool in proving the 
transcendence of e. Later, his student Pade systematically studied the scalar case p = 1 
[30] and such approximants are now called Pade approximants. 

From our perspective interpolating at infinity is more convenient than interpolating at 
the origin in the following sense. Any holomorphic function can be written as a Cauchy 
integral of its boundary values on any curve encircling a domain of analyticity. For a holo¬ 
morphic function on a domain we will use the terminology trace to mean the boundary 
values of the function on the boundary of the domain. When the function is holomorphic 
at infinity, such an integral representation can in some important cases be deformed an¬ 
alytically into an integral over a "one dimensional" set. A particular fruitful example of 
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this principle is when we deal with Markov functions. Let p? be a positive measure and 
take 


(i- 3 ) 


fi(z) 


dpdx) 

x —z 


supp(pi) C [at,biJ c 1R. 


The vector f is called an Angelesco system if [apbj] n [at,bt] = 0 for ) ^ i (such systems 
were initially considered by Angelesco [1] and later rediscovered by Nikishin [27]). Ange¬ 
lesco has shown that Qp has exactly nt zeros on [at, bj. This means that an Angelesco 
system is an example of a perfect system, i.e., a system for which every multi-index is 
normal, that is, deg(Qp) = | ft |. As far as the asymptotic behavior of 7ijJ' is concerned, 
convergence properties are not as straightforward as one could hope. Given an Angelesco 
system and a sequence of multi-indices such that ut/| ft t —t ct > 0, ci + • • • + Cp — 1, 
Gonchar and Rakhmanov [21] have shown that for each j the complement of U? , [at,bt] 
is separated by a system of analytic arcs into two domains, say DT, containing the point 

at infinity, and Dr, possible empty, such that rif converges to fj in Dt and diverges to 
infinity in D j”. Moreover, the polynomials Q n - can have an asymptotic zero distribution 
(in the sense of weak* convergence) on a strict subset of uL =1 supp(p.|) ( pushing effect). The 
pushing effect always implies existence of a divergence region but the reverse implication 
is not true. More detailed (strong) asymptotics for Hermite-Pade approximants to Angele¬ 
sco systems when p = 2 and the weights d iq/dx satisfy the so-called Szego condition, was 
obtained by the first author in [2]. 

Another class of Markov functions for which positive convergence results were ob¬ 
tained is now known as Nikishin systems [28]. The functions t? in (1.3) form a Nikishin 
system if they all are supported on the same interval [a, b] and the Radon-Nikodym deriva¬ 
tives dpj/dpi, j e { 2 ,.. .,p}, form a Nikishin system of order p — 1 with respect to some 
interval [c, d] such that [c, d] fl [a, b] = 0. Nikishin himself [28] has shown that such sys¬ 
tems are perfect when p = 2 and the Hermite-Pade approximants converge uniformly 
outside of the interval [a, b] in this case. This puts Nikishin systems more in line with the 
Pade case p = 1 (Markov theorem [26]) as neither the pushing effect nor the possibility of 

non-empty divergence regions appears for them. However, Nikishin systems do possess 

f 2) 

one new phenomenon, namely overinterpolation. It turns out that Re" has zeros on [c, d] 
that are dense on this interval. It took 30 years to prove that Nikishin systems are perfect 
for any p [17]. In [17], Fidalgo Prieto and Lopez Lagomasino also proved uniform conver¬ 
gence for multi-indices close to the diagonal. Strong asymptotics in the case of diagonal 
multi-indices and Szego weights was derived by the first author in [3]. 

It is interesting to observe that the first result on strong asymptotics of Hermite-Pade 
approximants was obtained by Kalaygin [24] for the case of two touching symmetric inter¬ 
vals (the limiting case of an Angelesco system). 

As often happens in mathematics, the treatment of Angelesco and Nikishin systems can 
be unified under the umbrella of generalized Nikishin systems (GN-systems) as introduced 
in [22], where Gonchar, Rakhmanov, and Sorokin defined a system of Markov functions 
with the help of a rooted tree graph and considered the question of uniform convergence 
of Hermite-Pade approximants to such a system. In such a set-up, an Angelesco system 
corresponds to a tree where the root is connected by p edges to p leaves and a Nikishin 
system corresponds to a tree in which every node except for the final leaf has exactly 
one child. Strong asymptotics of Hermite-Pade approximants to GN-systems of Markov 
functions generated by more general (than rooted tree) graphs (admitting cycles) was 
derived by the first author and Lysov in [8]. An example of GN-system from [8] is a pair 
of two Markov functions in (1.3) where the support of one of them is strictly included in 
the support of the other, i.e., supp(g2) = [ a 2/b2] C supp(m) = [ai,b-|] and the Radon- 
Nikodym derivative dp2/ / dpi along [a 1, b 1 ] is a Markov function with support la?,b3 1 
where [a3,b3] fl [ a 1, b 1 ] = 0. Weak asymptotics of Hermite-Pade approximants to this 
example of Markov functions was derived by Rakhmanov in [32]. 
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We emphasize that in all the results we listed above the geometry of the problem is real, 
i.e., the supports of the limiting distributions of the poles of the approximants and the 
overinterpolation points belong to R. 

In this paper we consider Hermite-Pade approximants to a pair of Cauchy transforms 
(Markov functions) of generally speaking complex measures with overlapping supports 
and aim at strong asymptotics. The set up does not fall into the framework of GN-systems, 
even when the measures are positive, as their supports overlap (for GN-systems the sup¬ 
ports are either disjoint or one coincides with or contains the other). It turns out that both 
phenomena, the pushing effect and overinterpolation, appear in this case. Hermite-Pade 
approximants to a pair of Markov functions with overlapping supports were first consid¬ 
ered by Stahl with the goal of proving weak asymptotics [34, 33, 9]. He had the important 
insight that the geometry of this problem is complex, i.e., the overinterpolation points are 
distributed on analytic arcs in the complex plane (later, a similar effect was observed in 
[5]). This discovery was very unusual at the time because the input geometry (i.e., the 
supports of the measures generating the Markov functions) is completely real. Unfortu¬ 
nately, Stahl's results have never been published. This work was strongly motivated by 
the desire to provide a detailed proof of his findings (in an even more delicate setting of 
strong asymptotics). 

Opting here for complex measures is natural from the point of view of complex analysis. 
However, many techniques, like those in [2, 3], do not apply as they use positivity in an 
essential way. An approach that does not rely on positivity was outlined by Nuttall in 
his seminal paper [29]. There Nuttall conjectured that the main term of the asymptotics 
of Hermite-Pade approximants is a function solving of a certain explicit boundary value 
problem on some unknown Riemann surface. He identified this surface only in a handful 
of special cases. Elaborating on Nuttall's approach, the first two authors and Kuijlaars [7] 
pinpointed the algebraic equation which defines the appropriate Riemann surface in the 
case of two Cauchy transforms of complex measures supported on two arcs joining pairs 
of branch points in the complex plane (the simplest example is a complexified Angelesco 
system) and derived formulae of strong asymptotics in the case when the Riemann surface 
has genus zero. Below, we build upon the ideas developed in [7] and extend the results of 
[7] to the cases when the appropriated Riemann surface has positive genus (elliptic and 
ultra- elliptic case). 

In Section 2 we identify the Riemann surface in Nuttall's program by an algebraic 
equation, discuss its realization as a ramified cover of C, and construct a certain function 
on this surface whose level lines will geometrically describe convergence and divergence 
domains of the approximants. In Section 3, we construct the Nuttall-Szego functions 
that will provide the leading term of the asymptotics of the Hermite-Pade approximants. 
Finally, in Section 4 we state the main result of this work. The remaining part of the paper 
is devoted to the proofs of all the stated results. 

2 RIEMANN SURFACE 

Let a £ ( 0 , 1 ) be given. Our goal is to investigate Hermite-Pade approximants to a pair of 
Markov-type functions generated by measures with supports [— 1 , a] and [—a, 1 ]. To this 
end we consider the algebraic equation 

(2.1) A(z)h 3 — 3 B 2 (z)h — 2 B-| (z) = 0 , 

where the polynomials A(z), B 2 (z), and B-| (z) are defined by 

( A(z) := (z 2 - 1 )(z 2 -a 2 ), 

< B 2 (z) := z 2 p 2 , 

B,(z) := z. 
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for some parameter p > 0 . Denote by h^, k £ { 0 , 1 , 2 }, the three distinct branches of the 
algebraic function h. determined by (2.1). Naturally, these branches satisfy 


(2.2) 


if ho + h^ + h .2 = 0, 

hohi + holy? + h] h.2 = —3B2/A, 
k h 0 h 1 h 2 = 2 Bt/A. 


Hence, it necessarily holds that h^foo) = 0 and therefore h^fz) = c^/zH— • as z —> 00. 
It readily follows from the above equations that the constants are the solutions of 
0 = c 3 — 3 c — 2 = (c + 1 ) 2 (c — 2 ). Thus, we put 

f h 0 (z) = 2 -\ - 

(2.3) < ? as z-)oo 

hi(z) =-1- 

v z 

for i £ { 1 , 2 }. It can easily be checked that all three solutions of (2.1) are real for positive 
large x. Hence, we can label the branches so that for all x > 0 large enough 


(2.4) ho(x) > hi (x) > h 2 (x). 

Denote by 93 the Riemann surface of h. It is a three-sheeted ramified cover of C. We 
shall denote by 93 ^, k £ {0,1,2}, the sheet on which h coincides with hj- (a particular 
realization of this surface is specified in Theorem 1 below). We denote by z a generic point 
on 93 with natural projection n(z) = z £ C. If we want to specify the sheet, we write z^ 
for z^ £ 93 ^ so that 7t(z^) = z. Thus, h^fz) = h(zf k ^). Consider the differential 


(2-5) 


dhf(z) := h(z)dz. 


where z is a generic point on 93 . The choice of the parameter p is driven by the following 
condition: 


(2.6) N(z) := Re 



is a well defined ( single-valued ) harmonic function on 93 . 


When it exists, N(z) is defined up to an additive constant. If we denote by the restric¬ 
tion of N to 93 then it is easy to see that Ni + N2 + N3 is a well defined harmonic 
function in C and therefore it is constant. Thus, we normalize N so that 


(2.7) 


Ni(z) + N 2 (z) + N 3 (z) = 0 , zee. 


Notice also that (2.6) is equivalent to requiring that all the periods of the differential dN 
are purely imaginary. 


Theorem 1 . Consider the algebraic equation (2.1) with a e ( 0 , 1). 

(I) If a e ( 0 , 1 /y/2), then there exists p £ (a, y/[\ + a 2 )/ 3 ) such that condition (2.6) is 
fidfilled. In this case 93 has eight ramification points whose projections are {± 1 , ±a} and 
{±b, ±ic}/or some uniquely determined b £ (a,p) and c > 0 . Moreover, the surface can 
be realized as on Figure i(a); 

(II) If a = l/y/ 2 , then condition (2.6) is fidfilled for p = 1 /y/ 2 . In this case 93 has four 
ramification points zvhose projections are { ± 1,±1 /y/ 2 .) and it can be realized as on 
Figure i(b); 

(III) If a £ [f/y/ 2 , l), then condition (2.6) is fulfilled for for p = y/{T+ a 2 )/ 3 . In this case 
93 has six ramification points zvhose projections are (± 1 , ±a} and {±b}/or some uniquely 
determined b £ (p, a). Moreover, the surface can be realized as on Figure i(c). 

The points b and c in Case I and the point b in Case III can be explicitly computed as 
they are solutions of a certain explicit quadratic or linear (in z 2 ) equation whose parame¬ 
ters depend on a and p. 

It follows from Theorem 1 that 93 has genus g = 2 when a £ ( 0,1 / y/ 2 ) and genus g = 1 
otherwise. Moreover, in Cases I and III, all the ramification points have order 2 while in 
Case II the points ±1 have order 2 while ±1 /y /2 have order 3. 
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9 i (0) o 9 t (2 > 9 i (0) o 9 i (1 > 9 i (0) o 9 i (2 > 


(b) Case II 


— 1 —a —b b 

9? 0) ^ JR*) 9l(°) ^ 91^ * 

(c) Case III 


a 1 

fR<°) ^ SK*) 


Figure l. Ramification points of 91 and the cuts (black curves) along which the 
sheets 91 ( 0 ) , 9 t m , 9 t (2) are glued to each other. 


In the light of Theorem l, it will be convenient to fix the following notation: 


(2.8) 


Ao 
Ai 
A 2 
A 2 i 
- A 22 


7t(cycle that separates 91 ( 11 and 91 * 2 '), 

7t(cycle that separates 91 ^ 0) and 91 * 1 '), 

A 2 i U A 22/ 

7t(the left cycle of the chain that separates 91 10 and 91 12: ), 
7t(the right cycle of the chain that separates 91 *°' and 91 f 2: . 


Clearly, Aq is defined only in Case I. 

Let dhf be defined by (2.5). The function 


(2-9) 


®(z):= exp 


d]\f 


is analytic, except for a double pole at oc/ 01 , and multiplicatively multi-valued on 91 . 
Moreover, it is single-valued in 91 a p := 91 \ Ut =1 ( a i U pj), where {at, is a homol¬ 

ogy basis on 91 . Later on, see Figure 2, we shall specify the basis in more detail, but right 
now it is sufficient to assume that each cycle y G {at, J 3 t }? =1 possesses a projective involu¬ 
tion: |y n 7t ' (z)| = 2 for any z £ 7t(y) which is not a branch point of 91 (the involution 
is then defined by mapping a point on y to the other one with the same projection). We 
normalize ® so that 


(2.10) 


(X)( 0) cD (1 )<D f2) = 1 in C, 


where < 1 * : k 1 is the pullback to C of the restriction of ® to 91 ' k '. Let us show that such 
a normalization is indeed possible. Since 91 a p is simply connected and dhf has integer 
residues, the restriction of ® to 91 a p is single-valued. It satisfies 

{ exp{ 27 tio)ij on a\, 
exp{27tiT|} on p|. 


1 ^ i ^ g, where the constants and X{ are real (this is guaranteed by (2.6)) and given 
by 


oh 


—— | dN, 
2m Jfh 


Ti ' 2 tii 


dW 


(2.12) 
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Pi p2 




(b) Cases II and III (there is no gap between the cuts on 9 t (0 ' in 
Case II, but this does not affect the choice of the cycles) 


Figure 2 . Homology basis for 91. 


At times, it will be convenient to use vector notation x = (xi , X2) T and cu = (up, up)" 1 " if 
g = 2 and x = (xi) and cu = (up) if g = 1 . Since O has a double pole at 00 and simple 
zeros at 00 ( 1 ), 00 i 2 i, we can write 

f O (0 i(z) = 


( 2 . 13 ) 


a>W 


C 0 z 2 + ■ 


= Ciz ” 1 +• 


as z —> 00. 


It follows from (2.7) that log |<Di 0 i<D* 1 )(])( 2 )| = 0 on C. Hence, if we choose the normaliza¬ 
tion Co Ci C 2 = 1, then (2.10) is fulfilled due to our choice of the homology basis: since 
the cycles possess projective involutions, exactly two pullbacks of <t> have jumps at each 
point belonging to the projection of a homology cycle, moreover, the jumps are reciprocal 
as the projected parts of the cycle coincide as sets but have opposite orientations. 

Due to the symmetries of the Riemann surface, the vectors cu and x take special forms. 
To be more specific, let us fix a homology basis. We choose 7r(( 3 1) = A21 and, in Case I, 
7t((3 2 ) = At, see (2.8), while the a-cycles are as on Figure 2(a) in Case I and the cci-cycle 
should be chosen as on Figure 2(b) in Cases II and III. 


Proposition 2 . Let {op, |3 i }® =1 be the homology basis which we just fixed. In Case 1 , one has 

(2.14) cu = (tu,2(l — cu)) T and t=(t, —■x) T 

for some real constants x = t(cl) and cu = cu(a) £ ( 1 / 2,1 ). In Cases II and III, one has 

(2.15) cu = x=(l/2). 
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Recall that |®(z)| = exp{N(z)} is a single-valued function on 91 , see (2.7). The asymp¬ 
totics of the Hermite-Pade approximants will depend on the relative sizes of the different 
branches of |<T>| (in other words, on the size of the branches of N). To this end, define Iltj k 
to be open subset of C such that 

O ijk := {z: |0> (i) (z)| > |0 (j) (z)| > |0 (,c) (z)|}. 

We also define the closed set F := Foi U Fq 2 U Fi2 by 

r tj :={z: |0 (i Hz)| = 

Clearly, C \ F = U iyj/kyi tiijk- Then the following theorem holds. 

Theorem 3. In Case I, the regions £lij k are distributed as on Figure 3(a); in Case II, the domains 
are distributed as on Figure 3(b); in Case III, the domains are distributed as on Figure 3 (c,d). 



(d) Case Illb 

Figure 3. Domains Dij k . 


We prove Theorems 1 and 3 as well as Proposition 2 in Section 6. 

3 NUTTALL-SZEGO FUNCTIONS 

To define the Nuttall-Szego functions, we first need to formulate a certain Jacobi inversion 
problem. To this end, denote by 

d(z) = (Q 1 (z),...,Q g (z)) T 










3 NUTTALL-SZEGO FUNCTIONS 


the vector of g normalized holomorphic integrals on 91 . That is, £l(z) is a vector of analytic 
and additively multi-valued functions on 91 that are single-valued on 91 and satisfy 

__ r, ns 

Be, on at, 

on 0t, 


(3-i) 


a+ - a - = 


B := 


O dQ; 


Pi 


j,i=1 


where eq is the i-th vector of the standard basis in ]R 9 . it j s known that B is a symmetric 
matrix with positive definite imaginary part when g = 2 and is a complex number with 
Im(B) > 0 when g = 1 . The Jacobi inversion problem we need consists of finding an 
integral divisor of order g, say I) n , such that 


(3.2) Q(B n ) = O(goo^) + c p + n(cu + Bt) (mod periods dO.), 


where c p is a vector of constants that depends only on the approximated functions and 
is defined further below in (7.11), t and d 3 are the vectors defined after (2.12), and the 
equivalence of two vectors c, e € C 9 is defined by 

(3.3) c = e (jnod periods dO^ <t=> c-e=j+Bm, j,rheZ 9 . 

It is known that the Jacobi inversion problem has a unique solution when g = 1 . Hence, 
Bn is well defined for all n in this case. Moreover, the following proposition holds. 


Proposition 4 . Assuming g = 1 , let B n be the unique solution of (3.2). Then I>2m = ®o 
and T>2m+i = Bi f or all m ^ 0 . Moreover, either Bo 7^ oo* 0 ^ or B] 00^ and therefore 
N* := {n G N : B n 7^ 00^} is infinite. 

The unique solvability of a Jacobi inversion problem is no longer guaranteed when 
g = 2 . However, it is known that solutions are either unique or given by any special 
divisor 1 2 . Notice that integral divisors of order 2 can be considered as elements of 91 2 /Il 2 , 
where IB is the symmetric group of two elements. This is a compact topological space. 
Hence, we can talk about limit points of sequences of integral divisors. 


Proposition 5 . Let g = 2, then (3.2) has a unique solution for at least one of the indices n — 1, n 
for any n e N. Moreover, there alzvays exists an infinite subsequence N* such that no limit point 
o/{B n } ng]N is a special divisor or is of the form 00 ^°^ + w for some w G 91. 

The following theorem is the main result of this section. 


Theorem 6 . Let pi and p2 be functions holomorphic and non-vanishing in a neighborhood of 
[—1,1]. In Case I, assume in addition that the ratio P2/P1 holomorphically extends to a non¬ 
vanishing function in a neighborhood of A$. Denote by A the chain that separates 9 f l °), 9 t (1 ), 
and 9 f^ 2 ' (n[A) is a union of sets defined in (2.8)]. For each nfor which (3.2) is uniquely solvable, 
there exists a function, say TV, meromorphic in 91 \ A such that the zero/pole multi-set ofW n is 
given by 

(3.4) B n -I- (n+ l)(oo^ ^ + 00^) — 2 noo^. 

Its traces are bounded except at the branch points of 91 where TV behaves like 

(3.5) |TV(z)| ~ |z-e| -1/4 as z -> e. 


1 Recall that an integral divisor of order d is a formal expression D = JV mzi, where rn £ IN and JV m = d. A 
principal divisor is an expression ^j c =1 ruti — £T =1 mi"Wt such that there exists a rational function on 9 t with 
a zero of multiplicity ru at ti for each i 6 {1 ,..., k), a pole of order mi at Wj for each i £ {1,. .., j }, and is 
otherwise non-vanishing and finite (in particular, one has that 2Zu=l n t = H[=i mi). 

2 On genus 2 surfaces, a special divisor is an integral divisor of order 2 such that there exists a rational function on 
9 t with simple poles at the elements of the divisor (a double pole if they coincide) and otherwise regular. 
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for every branch point e of IK except when e = ±1 /\fl (Case II only) where the exponent — 1/4 
should be replaced by —1 / 3 . Furthermore 




= ± 

(<>; 

f 

pi 

on 

A° 

^1' 

« 

f 

= =F 

(vr 

f 

P2 

on 

A° 

a 21 ' 

k 1 : 

f 

= ± 

(<>; 

f 

P2 

on 

A° 

^22' 

k 1 : 

f 

= ± 


f 

(P2/Pl) 

on 

> 

O ° 


where A° is the interior of the arc A. Moreover, ifW is afiinction meromorphic in 93 \ A satisfying 
(3.6), (3.5), and (3.4) with T) n replaced by any other integral divisor D, then M 7 is a constant 
multiple ofW n . 

It might seem, especially after looking at (3.2), that the construction of the Nuttall-Szegd 
functions M 7 tl , particularly, the divisors T) n , depends on the choice of the homology basis. 
However, this is precisely the point of the uniqueness part of Theorem 6 that the functions 
'F tl are independent of the intermediate steps of their construction. 

The Nuttall-Szego functions T tl are certain multiples of the functions <l> n . In fact, their 
ratios form a normal family on each 93 £ obtained from 93 by excising circular neighbor¬ 
hoods of radius £ around the branch points. 


Theorem 7. Let N* C N be a subsequence as in Propositions 4 or 5. For each z > 0 , there exists 
a constant C £ (N*) > 1 such that 


|Vnl<C e (N*)|® n | in 93 £ , 

|¥ n | /i CaN*)" 1 |O n | in 93^ n tt -1 {\z\ > 1/e}. 


We prove Propositions 4 and 3 as well as Theorems 6 and 7 in Section 7. 


4 ASYMPTOTICS OF HERMITE-PADE APPROXIMANTS 


Below, we consider vector functions f := (f 1, If) of the form 


( 4 -i) 


fj(z) 


1 

27 ti 


P)( 


p. x —z 

1 ^ 


dx, j G {1,2}, 


where Fi = [—1, a] and F2 = [—a, 1], a G (0,1), and pj are holomorphic and non-vanishing 
in a neighborhood of [— 1 , 1 ]. Additionally, we impose the following condition on the 
functions py 


Condition 8 . The ratio P2/P1 extends from (—a, a) to a holomorphic and non-vanishing function 

• in a domain that contains in its interior the closure of all the bounded components of the 
regions 12 ^^ in Case I, see Figure 3(a); 

• in a domain whose complement is compact and belongs to the right-hand component of 
C>021 in Cases II and Ilia, see Figures 3 (b,c); 

• in the extended complex plane, i.e., the ratio is a non-zero constant, in Case Illb. 


Condition 8 needs to be regarded in the following context. As conjectured in [4] and 
supported by the results for Nikishin systems, given f of the form (4.1), the appropriate 
Riemann surface 93 must depend on the analytic continuation of P2/P1 from F( n F2. In 
this work, on the other hand, we fixed the surface by considering (2.1), which necessitates 
a condition on the continuation of P2/P1 of the above type. 

In what follows, we assume that a € ( 0 , 1 ) is fixed and the vector f = (f-|, f2) is given 
by (4.1), where pi and P2 satisfy Condition 8; 'l ; Tl are the Nuttall-Szego functions of 
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Theorem 6 corresponding to a and the functions pi and p2, constructed on the Riemann 
surface 9 F of h realized as in Theorem 1; the constants C n are defined by 

(4-2) C" 1 := lim ^ 0) (z)z- 2n ; 

z —too 

Qft are the monic denominators of the Hermite-Pade approximants tt/i to f, which are 
also multiple orthogonal polynomials with respect to the weights pi and p2, see (5.1)), 
and R ^ 1 , i £ { 1 , 2 }, are the linearized errors of approximation corresponding to the multi¬ 
indices it = (u,n), n £ N, see (i.i)-(i.2). 

To describe the asymptotics of r| (', it will be convenient to put 

H := rf 2 \A 0 , [ Hi := (^021 un 2 oi) \ (Ti2 U A] U A 2 ), 

, and < _, _ _ 

r 2 := r+, (n 2 := (q £ 21 uaj 01 ua) \ (r+ ua, ua 2 ), 

where Ai and A2 were introduced in (2.8), F^ := Pi2 n{±Re(z) > ^ijk n 

{±Re(z) > 0 }, and £1 is present only in Case I and is equal to the unbounded component 
of Qq-\2' see Fig ures 4 and 3. Observe that II1 = CI2 = 0 in Case Illb. Define 

4i 2) in II 1 , \ ^ 2) in0 2 , 

and Wn := < 

1 otherwise, [ pf^^ 1 otherwise. 

In Case I, it follows from (3.6) that 1 is the analytic continuation of across Ao until 
Fj_, which always exists by the analyticity of pi's. In Cases II and III it still holds by (3.6) 
that is an analytic continuation of M'j} . However, this time one needs to continue T/ 1 1 
through A^ into 91 ■ 01 and then through Af) into 91■ 2 while M'/ - needs to be continued 
through A 22 into 91 f 0: and then through A^ into 91 11 1 . 

For any 5 > 0 , we further define 

D\Tg :={z: |lm(z)| < 6, Re(z) £ Ai UA2, dist (Re(z),{endpoints of Ai and A2}) > 6}, 


( 4 - 3 ) 


:= 


.£1 

P2 


{ {z = t + x : xe (—6,6), t £ Fi, 6 < |lm(t)| < c — 6} in Case I, 

{z = t + x : x £ (—6,6), t £ Fi, 6 < ]lm(t)| < 1 / 6 } in Case II, 

{z : dist(z, Ti) < 6, |lm(z)| > 6 when |Re(z)| < 1 } in Case III. 


Finally, let us introduce the following notation. Given a sequence of functions F n , a 
sequence of finite multi-sets X n , and positive numbers e n , we write 


Fn — 0 (e n ;X T 




• n 

xGX n 


z-x 


VTl + W 2 )0 + M ; 


= 0(e T 


Then the following theorem holds. 

Theorem 9 . Let N* be a subsequence from either Proposition 4 or Proposition 5 (depending on 
whether a ^ 1 /y /2 or a < 1 and i £ { 1 , 2 }. Then 


( 4 - 4 ) 


Qn = CnWf 1 (1 +0(en;Xn)), 
r £ 5 - Cn^ i] (l-FO(en;X^)), 


locally uniformly in C \ (Ai U A2) and C \ (Fi U Ft), respectively, where X n is the multi-set of 

zeros of^n * in C \ (Ai U A2), X^' is the multi-set of zeros ofw[^ ] in C \ (Ft U Ft), and 3 e n —t 0 . 

Moreover, there exists 6q > 0 such that for all 6 £ ( 0 , 6q) one has 

( 4 - 5 ) 


Qrt = 


Cn<°i 0 + 0(en;X n+ )) + C n <°! (1 + 0(en;X n _)), 




= C 


(l +0 (e n ;X£]_)) +Cn^l (1 4-0 (e n ;X£I)), 


n € N*, 


3 In Cases I and III one has e n = n 1 and e n = n 1 1 6 in Case II. 
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Figure 4. The sets (non-horizontal bold lines), the domains 0 ; (shaded re¬ 
gions), and the boundary of the domain of holomorphy of ' (bold lines). 


locally uniformly in N5 and Ng 1 *, respectively, where are functions holomorphic in N5 that 
coincide with 0 n A'\ U A2 and X n ± are their multi-sets of zeros in Dvfg, and are 

functions holomorphic in \ F, that coincide with ^^ on F° \ F t and xj^j. are their multi¬ 
sets of zeros in Xlg 1 *. 










5 MULTIPLE ORTHOGONAL POLYNOMIALS 


It follows from Theorems 3 and 7 that in Cases I and III we observe the pushing effect, 
i.e, a b, and the presence of divergence regions (both phenomena are observed in 
Angelesco systems). Indeed, according to (1.2) and (4.4), one has locally uniformly in 

c\([-l, 1]urt) that 

f _ (i) _ 1 +0(en;Xn 1 ) 

' 1 + 0 (e n ;X n ) ‘ 

Hence, the error of approximation C — 7t/ 1 is geometrically small in (|z| ^ 1 /e} for all £ 
small enough. On the other hand, in Case I, the error f 1 — 7tjJ 1 is at least geometrically 
big on compact subsets of £ 1 J 01 , see Figure 3(a), (it could be infinite if the elements 
of D n belong to 10 and project onto this component) and the error f2 — 7tL - ' is at 
least geometrically big on compact subsets of £l/ 01 ■ In Case III, both components of the 
approximant diverge in O-|02 and in Case II there are no divergence domains. 

As to the zeroes of the functions, it can be deduced from Theorem 3 and (4.5) that Qfj 
must vanish in 3 \f§ and R/ must vanish in hi/ 1 *, which is precisely the phenomenon of 
overinterpolation first observed in Nikishin systems. 


5 MULTIPLE ORTHOGONAL POLYNOMIALS 


The basis of our approach to asymptotics of Hermite-Pade approximants lies in their 
connection with multiple orthogonal polynomials. It is quite simple to verify that if the 
functions ft are of the form (4.1), then (1.2) is fulfilled if and only if 


(5-i) 


Q rt (x)x k pi(x) dx = 0 , k e { 0 ,.. .,ut - 1 }. 
U 


Moreover, the linearized error functions R 

,(ii. . 1 


(t) 


(5- 2 ) 


Rl 


z = 


2 m 


U 


admit the following integral representation: 

QrtMpi(x) 
x —z 


The analysis of the system (5.1 )—(5.2) then proceeds via its reformulation as a matrix 
Riemann-Hilbert problem. This fundamental fact in the theory of orthogonal polynomials 
was first revealed by Fokas, Its, and Kitaev [18, 19] and the extension to multiple orthogo¬ 
nal polynomials was given in [37]. Set hi := (n— l,n) and fi2 := (n,n— 1 ), and assume 
that the index n is such that 


(5-3) 




+ ■ 


for some constants m/C Under condition (5.3), the matrix 


( 


(5-4) 


Y := 


’ Qftj 


4 1] 

(Un(U 


mV ’R 


n-1 


m (2) n m (2) p (l 

\ Ln n Q112 Tiin R 


H-2 


Rft ’ 

"AX 1 


solves the following Riemann-Hilbert Problem (RHP-Y): 


(a) Y is analytic in C \ [— 1 , 1 ] and 

lim Y(z) diag ( z" 2n , z n , z 11 ) = I, 

z—>00 V / 

where diag(-, -, •) is the diagonal matrix and I is the identity matrix; 

(b) Yhas continuous traces on (— 1 , 1 ) \{±a} that satisfy Y + = Y_J(xi Pi ,XzPz), where 

1 x y \ 

0 10 
0 0 1 / 

and Xi is the indicator function of Ft; 
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(c) Y(z) = O flog |z — e|) as [— 1 , 1 ] ^ z —> e for e e {± 1 , ±a} 4 * . 

Vice versa, if RHP-Y is solvable, then the solution is necessarily of the form (5.4) and 
(5.3) holds. To prove Theorem 9, we then follow the framework of the non-linear steepest 
descent method for matrix Riemann-Hilbert problems, first introduced in the 2 x 2 case 
by Deift and Zhou [14]. The proof of Theorem 9 is carried out in Sections 8-10. 


6 GEOMETRY 

This section is devoted to proving Theorems 1 and 3 together with Proposition 2. In 
Section 6.1 we establish that 91 can indeed be realized as in Figure 1. In Section 6.2 we 
justify the the choice of the parameter p, thus, finishing the proof of Theorem 1, while 
simultaneously proving Proposition 2. Finally, we prove Theorem 3 in Section 6.3. 


6.1 Realization of 9 \ 


It can be readily computed that the discriminant of (2.1) is equal to 
D(z) = 108 A(z)[B|(z)-A(z)Bf(z)] 

= 108 A(z) [(1 + a 2 — 3 p 2 )z 4 + ( 3 p 4 — a 2 )z 2 — p 6 ]. 
Assume first that we are in Case I, i.e.. 


(6.1) 


(6.2) 


a G 0, 


1 

71 


and p £ a. 


1 +a 2 


Since the polynomial B| — AB 2 , which is symmetric and of degree 4, is negative at the 
origin and has positive leading coefficient, it follows that it has four zeros, which we 
denote by ±b and ±ic, where b, c > 0 . Furthermore, be (a,p) since B| — ABy is positive 
at p and negative at a. Observe that if a point is a branch point of h of order 3 (all three 
branches coincide) and h is finite at this point, then necessarily all three branches are equal 
to zero there. Hence, ±b, ±ic are branch points of order 2. Furthermore, the first equation 
in (2.2) implies that neither of the points ±l,±a can be a pole of one the branches while 
the second equation implies that the branches cannot have a cubic root singularity there. 
Hence, the points ± 1 , ±a are branch points of order 2, two branches are infinite at them 
and one is finite. 

It follows from the above discussion that we can analytically continue the branches h k 
so that the inequalities in (2.4) hold for x > 1 . At 1 two branches blow up and, of course, 
all three add up to zero. This is possible only if ho tends to 00, h 2 tends to —00, and hi 
remains bounded. Thus, 1 is a branch point joining 93 * and 91 21 . It can easily be seen 
from (2.1) that all the branches must satisfy 


( 6 - 3 ) 


Mz) = h k (z). 


Therefore, all the branch cuts must be conjugate-symmetric. Thus, the branch cut starting 
at 1 must end at b. That is, ho± = h 2 =i= on (b, 1 ), which immediately implies 

(6.4) (h 0 -h 2 )+ + (h 0 -h 2 )_ = 0 


on (b, 1 ). Moreover, (6.3) implies that the traces above are purely imaginary. As ho, h 2 are 
unbounded near 1 and bounded near b, one has 

(6.5) (h 0 -h 2 ) 2 (x) = ’j’fi(x), 

where fi is holomorphic, non-vanishing, and real on (a,00). As tq must be positive for 
x > 1 and is non-vanishing, it is, in fact, positive for x > a. Therefore, (ho — h 2 )(x) > 0 , 


4 In fact, in each case the entries of at least two columns remain bounded. However, the above simplification does 

not affect the forthcoming analysis. 
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x G (a, b) U ( 1 , oo). Moreover, since hi is real and non-vanishing for x > a and is negative 
for x large enough by (2.3), hi is negative for all x > a. Thus, 

(6.6) h 0 (x) > h 2 (x) > hi (x), xG(a,b). 

As in the case of 1, two branches are infinite at a and one is finite. The inequalities in 

(6.6) imply that the unbounded branches are ho and hi and therefore a is a branch point 
between 91 ^ and 9t 1 1 . The branch cut is (-Q, a) since it must be along the real axis. That 
is, ho± = hi =|= on (—a, a) and 

(6.7) (h 0 -hi)+ + (h 0 -hi)_ = 0 

on (—a, a) where the traces are purely imaginary. As ho, hi are unbounded near ±a, one 
has 

(6.8) (h 0 -hi) 2 (x)= 2° ^*2 ' 

where f a is holomorphic, non-vanishing, and real on (—b, b). The function f a is positive 
on (a,b) according to (6.6) and therefore is positive on (—b,b). Hence, (ho — hi )(x) < 0 
for x G (—b, —a). Since both branches blow up at —a, one has 

(6.9) hi (x) > h 2 (x) > ho(x), x G (—b, —a). 

Furthermore, it is quite simple to deduce from (2.1) that the branches h k must satisfy 
hk(— x) = —hj k (x) for all k G { 0 , 1 , 2 }, where jic G { 0 , 1 , 2 } depends on k. According to 
(6.6) and (6.9), the considered continuations of h^ satisfy = k for |x| G (a, b). Therefore, 
—b is a branch point for ho and h 2 as well as — 1 . In particular, (6.4) remains valid for 
x G (— 1 ,— b). Hence, 

(6.10) (h 0 — h 2 ) 2 (x) = —y-j-f 2 1 (x), 

where f_i is holomorphic, non-vanishing, and real on (d, —a) for some d < — 1 . From 
(6.9), we know that f_i is negative on (—b, —a) and therefore it is negative on (d, —a). 
Thus, (ho — h 2 )(x) < 0 for x G (d, — 1 ). This and the blowing up of ho and h 2 at — 1 imply 
that 

(6.11) h 2 (x) > hi (x) > h 0 (x), x G (d, — 1 ). 

On the other hand, write 

(6.12) ht(z) = —- + H- 

Z Z z 

for i G { 1 , 2 }. Plugging this expansion into (2.1) and considering the 1 /z term on the 
left-hand side, we get that 

-> 1 + a 2 -> 

(6.13) 0Cr = — - p 2 >0. 

This means that the inequality hi (x) > h 2 (x) holds for all |x| large enough. That is, 

hi (x) > h 2 (x) > h 0 (x), x G (-00, d'), 

where necessarily d' < d. As there are no branch points between d' and d, there should 
be a branch cut passing between them and this cut should necessarily be between 91 ^' 
and 9 t^ 2 '. In other words, ±ic are branch points of hi and h 2 . This finishes the proof of 
the claim that 91 can be realized as in Figure i(a) in Case I. 

Assume now that we are in Case II, i.e., 

Q = P = 1 /VI. 

Then the discriminant of (2.1) is equal to D(z) = 27 (z 2 — l)(z 2 — 1 / 2 ). As before, ±1 are 
branch points of order 2. Furthermore, we know from the third identity in (2.2) that some 
of the branches are unbounded near ±l/\/ 2 . However, if the branching were of order 2, 
the left-hand side of the second relation in (2.2) would be unbounded near ±1 /V2, but it 
is bounded there. Hence, ±1 /y/ 2 . are branch points of order 3. 
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Plugging a power expansion for tp, i £ { 1 , 2 }, into (2.1), one can compute that 

(6.14) ht(z) =-f H-, 

z z A 

where the ( 3 |'s are solutions of ( 3 2 + (3 + 1/6 = 0 . In particular, (hi — h 2 )(z)z 3 = l /\/3 + 
0(1 /z] for all z large. Then, by repeating the initial steps of the analysis for Case I, we see 
that (2.4) extends for all x > 1 and that 


(6.15) h 2 (x) > hi (x) > h 0 (x), x€ (-oo,-1). 

Therefore, (— 1 ,- 1 / y/l) and ( 1 / y/l, 1 ) are the branch cuts for ho and h 2 . Around the point 
1 /y/l it can be directly verified that (2.1) is solved by the following Puiseux series: 

' / 1 N-I/3 / -i \ 1/3 

^( z ~7l) hV3(z]+£ ’ 2 ( z_ 7l) H V3(z) ' 

where £, is any solution of £, 3 = 1 , H(z) := (l + 2 zvT^z?)/( 2 z + y/l), and all the roots 
are principal. Since hi is negative and holomorphic for x > 1/V2, one finds that hi (z) = 
h(z; 1 ) locally around 1 /x/ 2 . On the other hand, 

ru u \ 1 < v f2 M 
(h 0 -h 2 ) (x) = —-, 

x z — 1 


(6.16) h(z;f.) = 




where f(x) is non-vanishing and holomorphic for |x| > 1 /y/l. Since f(x) > 0 for x > 1 , we 
get that f (x) > 0 for x > 1 / y/l and 

(6.17) h 0± (x) = h 2=F (x) = —^hi (x) Ti^===, xe(l/V / 2 ,l). 

Thus, ho+ has values in the fourth quadrant and ho_ has values in the first quadrant. As 
the first summand in (6.16) is dominant around 1 /y/l, we can conclude that 

(6.18) 

u , , f h(z;e 2 ™/ 3 ), lm(z) > 0, , w , f h(z;e 4 ™/ 3 ), lm(z) > 0, 

( h(z;e 4 ™/ 3 ), lm(z) < 0, \h(z;e 2 ™/ 3 ), lm(z) < 0, 

locally around 1 /yj 2 . From this it is easy to see that h 2 is holomorphic across the interval 
( — 1 / y/l, 1 / y/2) and this interval is the branch cut for hi and ho- 
Finally, assume that we are in Case III, that is. 


a £ 



and p 


1 +a 2 
3 


< a. 


In this case B} — AB 2 is a polynomial of degree 2 which has two roots ±b satisfying 
b £ (p, a). Exactly as in Case I, we see that all the branch points, namely {± 1 , ±a, ±b}, are 
of order 2. Furthermore, the same reasoning as in Case I gives that (6.4) holds on (a, 1 ). 
Since ho and h 2 are unbounded at both 1 and a, we get that 


(6.19) 


(h 0 -h 2 ) 2 (x) 


(x — 1 )(x — a) 7 


where fi is holomorphic, non-vanishing, and real on (b, 00). Since fi is positive for x > 1 
and hence for x > b, we can conclude that (ho — h 2 )(x) < 0 for x £ (b, a). The blowing up 
of ho and h 2 at a implies that 


(6.20) 


h 2 (x) > hi (x) > h 0 (x), x £ (b, a). 


It further follows from the third equation in (2.2) that two branches of h are negative and 
one branch is positive on (b, a). The inequalities in (6.20) show that the negative branches 
are ho and hi. This, in turn, implies that b is a branch point of ho and hi and so is —b. 
Hence, (6.7) holds on (—b, b). As all the branches are bounded at ±b, we have that 

(6.21) (h 0 — hi ) 2 (x) = (x 2 — b 2 )f^(x) 
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for x £ (—a, a), where fb is holomorphic, non-vanishing, and is negative on (b, a). There¬ 
fore, it is negative on the whole interval (—a, a) and we get that (ho — hi )(x) > 0 for 
x € (—a, —b). As before, the third equation in (2.2) implies that two branches of h are 
positive and one branch is negative on (—a, — b). Moreover, it also implies that h. 2 ( 0 ) = 0 . 
As h 2 (x) > 0 for x > 0 and it has no other zeros, it must be the negative branch. Thus, 

(6.22) h 0 (x) > h] (x) > h 2 (x), x € (—a, —b). 

Now, as the branches that meet at —a are unbounded, they must be ho and h 2 by (6.22). 
Hence, (6.4) holds on (— 1 ,— a) and we have that 

7 f-i M 

(6-23) (h 0 — h 2 ) 2 (x) = , 

(x+ 1 )(x + a) 

where f_i is holomorphic, non-vanishing, and positive on (—a, — b). As before, this means 
that f_i is positive for all x < —b and (6.15) holds in this case as well. This finishes the 
proof of the claim that 93 can be realized as in Figure i(c) in Case III. 


6.2 Choice of the Parameter 


We start with Cases II and III as they are much simpler. Here we show that for the choice 
of the parameter p as in Theorem 1 the condition (2.6) is fulfilled. That is, the period 
of the Nuttall differential dhf over any given chain on 93 is purely imaginary. The latter 
simplifies to showing that the periods of dhf on the cycles of a homology basis are purely 
imaginary. In fact, (2.15) implies that both periods of dN are equal to 7ti. 

In proving (2.15) we shall rely on the following observation: one has that h^f— z) — 
—hj k (z) for some j ^ C { 0 , 1 , 2 }, which can be deduced immediately from (2.1). In fact, we 
see from (2.4) and (6.15) that h^}— z) = — h^z) for all k € { 0 , 1 , 2 }. This implies that 


op = 


1 

27 ti 


A21 


(h 0+ — h 0 _)(x) dx = - 


2ni 


A 2 2 


(h 0+ -h 0 _)(x) dx. 


The sum of the last two integrals is equal to 1 by Cauchy's residue theorem applied to 
h 2 (recall that ho± = h 2 q; on A and h 2 (z) = — 1 /z + 0 (l/z 2 )). This gives the desired 
conclusion about cu 1. Furthermore, from the choice of the a-cycle, see Figure 2(b), one 
has for Case III that 


U = 


27 ti 


h 0 (x) dx + 


—b 


h 0+ (x) dx + 


—b 


h 0 (x) dx + 


h 2 (x) dx 


As ho is an odd function, the sum of the first and the third integrals is equal to zero. Since 
h 2 is odd as well, the fourth integral is zero too. Using the symmetry considerations once 
more, we can get that 

—b r b 


ho+(x)dx= ho_(x)dx. 

. b J—b 

Then by applying Cauchy's residue theorem to hi (recall that hi (z) — —l/z+ 0 (l/z 2 )), 
we see that 


U = 

In Case II, we have that 
U = 


1 

4 m 


-b 


J_ / r-b 
4 m 


h 0+ (x) dx + 


-b 


ho_(x) dx 


hi_(x) dx + 


-b 


hi + (x)dx = 1/2. 


1 

2ra 


-1/V2 

I /\/2 


h 0+ (x) dx + 


■1/V2 

-1/C2 


. h 2 (x) dx , 


and the conclusion r-\ = 1/2 follows from an analogous symmetry argument. This finishes 
the proof for Cases II and III of Theorem 1. 
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In Case I we again start with (2.14). The symmetry of the surface implies in this case 
that ho is an odd function, i.e., ho (— z) = —ho (z), the functions hi, i € {1,2}, are odd within 
the bounded domain delimited by Aq and Aq := {z : —z £ Aq}, and hi (—z) = — h.2(z) for 
z within the unbounded domain delimited by Ao and -Aq. Then it follows from (2.3) and 
the Cauchy residue theorem that 


2 = - 


1 

2m 

1 

2™ 

1 

271 


A 2 UAi 


(h 0+ — ho_)(x) dx 


a 2 


(h 0+ -h 2+ )(x) dx - — 


a 21 


(~f-i M) 


-b-x 

x+1 


dx + 


A, 

1 

271 


(It-oh 


-h 1+ )(x) dx 

faW 


Vz 2 —' 


= dx+ — 

2 271 


A 22 


fl M 


x-b 
1 -x 


dx. 


where we used the notation from (6.5), (6.8), and (6.10). The analysis right after (6.5), (6.8), 
and (6.10) implies that each of the last three integrals is positive. From the definition of 
our homology basis, it can easily be seen that the first integral is equal to oq, the second 
one is equal to o>2, and the third is equal to the first as ho is an odd function. That is, 
cut > 0 an d 2 o> j + m2 = 2 . The latter clearly implies the first relation in (2.14) and the fact 
that cu = cu-| £ ( 0 , 1 ). Applying now the Cauchy residue theorem to h.2, we get that 


1 = ^ 


1 

27 ti 


a 2 


2aq + -*-t 
2 m 


(h 2 + - h 2 _) (x) dx + — 
(h 2 + -h 1 + )(x) dx. 


(h 2 + — h 2 _)(x) dx 


where Aq is oriented from —ic to ic. Hence, we can conclude that aq £ ( 1 / 2 , 1 ), if we show 


that the last integral is negative. Since h.2 + 
can deform the path of integration to get 


h-| + is a trace of a holomorphic function, we 


1 

27 ti 


(h 2+ -h 1 + )(x)dx = - 


(6.24) 



+ 


+ 


(h 2 (t)-h 1 (t))dt 


(h.2(ix) — hi (ix)) dx 


(h 2 (ix) — hi (ix)) dx. 


where we used the symmetry h.2(— z) = —hi (z) to get the last equality. Notice also that 
the branches satisfy fq/z) = h-i/z) as follows directly from (2.1). Hence, 


h2 (ix) — hi (ix) = h.2 (—ix) — hi (—ix) = h.2 ( — ix) — hi ( — ix) = h.2 (ix) — hi (ix) 

for x > c. Hence, this difference is real there. Moreover, h_2(ix) — hi (ix) > 0 for all x large 
as follows from (6.12) and (6.13). Since the difference h2 — hi can be equal to zero only at 
the branch points of h.2 (ix) — hi (ix) > 0 for all x > c, which shows that the integral in 
(6.24) is negative as desired. 

To prove the second relation in (2.14), observe that our choice of the homology basis 
can be made so that the a-cycles are contained within the bounded domain delimited by 
Aq and — Aq, see Figure 2(a). Moreover, they can freely be deformed within the domain 
of holomorphy of dhf. Thus, it follows from (2.5) and (2.12) that 


27tiTi = 


—b 


ho (x) dx + 


— a 
•b 


h.2(x) dx + 


-b 


h.2 (x) dx + 


hi (x) dx + 


ho_(x) dx 


h 0 (x) dx + 
—27tiT2, 


h-2(x) dx + 


h-2(x) dx + 


hi (x) dx + 


h 0+ (x) dx 


where we used the fact hi/— z) = —h/z) with the bounded domain delimited by Ao and 
—Aq, which, in particular, implies that hi_(—x) = hi + (x) for x £ (—a, a). This finishes 




































6.2 CHOICE OF THE PARAMETER 


the proof of (2.14). It only remains to show that there exists a choice of the parameter p in 
(2.1) so that t from (2.14) is real. 

Showing that there is a choice of the parameter p £ (a, \/(l + a 2 ) / 3 j such that t in 
(2.14) is real, is equivalent to proving that 1 = 0, where 


(6.25) 


I := Re o h(z) dz 


«2 


In order to prove (6.25) we consider the limiting values of I = I(p) as p —> ^/(l + a 2 )/3 
and p —> a. If they have opposite signs, then clearly such a choice of p is indeed possible. 
We start with the case p —t •*/(! + a 2 )/ 3 . Define 


N k (z;p) 
Ni (z;p) 


Re (Jb M dx), h £ {0,2}, 
Re (Ja h i M dx ) + N o(a;p), 


z £ {Re(z), lm(z) ^ 0 }. 


Then it obviously holds that I(p) = N1 (ic;p) — N2(ic;p). In Section 6.3 further below we 
shall argue that 


(6.26) 


N(z;p*):= Re 


^ h.(x;p*)dx^, p* := ^/(l + a 2 )/ 3 . 


is a well defined harmonic function on the Riemann surface of h( ;p*) and that 

(6.27) N(oo^';p*) — N(oo^;p*) < 0 . 

As I(p) depends continuously on the parameter p, we can conclude that I(p*) < 0 . 
Let now p = a. In this case equation (2.1) becomes 

(6.28) (z 2 — 1 ](z 2 — a 2 ) h 3 — 3 (z 2 — a 2 ) h — 2 z = 0 , a 2 < 1 / 2 . 


The branch points with projections a and b of the curve (6.28) merge together into a triple 
branch point with projection a (this can be observed directly from (6.1)). To verify that 
1(a) > 0, we deform the cycle a2 = 0C2 (a) into a cycle a which is involution-symmetric 
and whose projection from say a, is as on Figure 5. That is, a emanates from a 


z(t_: 


1C 


•1 


a 


-a 0 a 1 

Figure 5. Contour cc for the limiting case p = a < I/a/2. 


into the first quadrant along some special arc z(t), t £ [0, t_], and then proceeds along the 
imaginary axis from z(t_) down to ic. 

As in (6.24) one has tq(z; a) = tu(z; a) and hi (—z) = — h.2(z) for all z large (and hence 
for all z £ (ic,ioo)), which leads to the same conclusion that (hi — h.2)(z) is real for 
z £ (ic,ioo) and therefore the real part of the integral of this difference on [ic,ioo) is 
equal to zero. Hence, we get that 


1 (a) = Re 


(hi (z; a) — h 2 (z; a)) dz 


Therefore, in order to complete the proof of Theorem 1, it remains to show that there exists 
an arc z(t), such that 


(6.29) 


z( 0 ) = a, Re(z(t_)) = 0 , lm(z(t_)) > c, 

Re((hi -h 2 )(z(t))z'(t)) > 0 , t £ ( 0 ,t_). 
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To construct this arc we use a parametrization of the algebraic curve (6.28). This parametriza- 
tion was suggested in [11] and it has the form 


(6.30) 

where 

(6.31) 


h = —i { £t + — ) C(t) 


^ = 1 , tec. 


Z = 


1-t 6 


(a_(t) +ia+(t)) 


a±(t) := ^±(1 - a 2 ) =F a 2 t ±6 , C(t) := 


a+(t) + it 6 a_(t) 
(1 -a 2 )(1 +t 6 ) ' 


Since the expression for h in (6.30) has the Cardano form 

h = f, A + f,” 1 B, A = —itC, B = —- C, 

the equivalence of (6.30) and (6.28) can be written (using the Vieta relations) as 

1 


-AB = C z = — 


z 2 — 1 

A 3 + B 3 =iC 3 (t 3 +t“ 3 ) = 


2z 


(z 2 — 1) (z 2 — a 2 ) 


With some work one can verify that plugging in the expression for C(t) from (6.31) and 
the expression for z(t) from (6.30) into the above equations produces an identity, which 
proves that (6.30) is indeed a parametrization of (6.28), (for details see [11]). Observe that 

(6.32) a±(t) VO, t e [0,t_], t_ := 



and that the parametrization (6.30) defines a Jordan arc with end points 

. ia , ia 2 

z 0 = a and z(t_ = , :, while ic = . ■ , 

VI - 2d 2 VI -2u 2 

As d > d 2 , the first line of (6.29) is satisfied. The local analysis around the point d 
shows that hi is given by (6.30) with £, = (i \/3 — l)/2 and for h 2 one needs to select 
£, = (—iV 3 — 1 )/ 2 . Then, after a tedious computation, we get 


(hi -h 2 )(z(t))z'(t) 


3V3t /1-2d 2 _.\ 

t 4 +1 2 + 1 \ u+ d_ / ‘ 


Hence, the second line of (6.29) follows from (6.32) as u 2 < 1 / 2 . This finishes the proof 
of Theorem 1 granted we can prove that N(z;p*) in (6.26) is a well defined harmonic 
function, which we do at the end of the upcoming subsection. 


6.3 The Regions n t j k 

We start with two general observations that are consequences of the single-valuedness of 
N(z) on 93 . First, the regions T>p k could be equivalently defined by 

O ijk := {z:Nj(z) >Ni(z) >N k (z)}, 

and this definition does not depend on the initial point of integration chosen in (2.6) as 
changing the initial point results in adding the same constant to all N k simultaneously. 
Secondly, let H be the analytic continuation of hi — h 2 from the point at infinity. Then H 
is an algebraic function and Re(H) is a well defined harmonic function on the Riemann 
surface of H. Hence, the set V is a subset of a projection of the zero level line of Re(H) to 
C. As such it cannot be dense in an open set. 

We start with Case I. It will be convenient to consider a slightly different realization of 
93 , namely, we shall suppose that 

(6.33) A0 = [—ioo, —ic] U [ic, ioo]. 
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Our first goal is to show that E in this case has the form as shown in Figure 6. It follows 



Figure 6. Case I: Pi 2 - solid lines, Pq 2 - dashed lines, and Pqi - dotted dashed line. 


from our first observation that lines Pq emanating from a projection of a branch point e 
of th can be described by 


0 = Re 



hj) (x) dx I , 


for z € Tij locally around e. That is, Pq is a trajectory of a quadratic differential jly — 
hj) 2 (x)dx 2 . The local behavior of the trajectories is well known, [35]. This implies that 
exactly one line of the set P emanates from the points ±1 and ±a ("hard edges") and 
exactly three lines (with angle n /3 between them) emanate from the points ±b, ±ic ("soft 
edges"). 

In order to proceed, recall the inequalities (2.4) and (6.6) as well as the decompositions 
(6.5) and (6.8). Furthermore, it follows from (6.3) and the fact that ho(— z) = —ho(z) that 
ho is purely imaginary on ilR. Analogously, we can conclude that Pq, i £ { 1 , 2 }, are purely 
imaginary on (—ic, 0 ) U ( 0 ,ic). 


0 

ho 


r 


Reh-o < 0 
Imh-o < 0 

a 

b 

Re h-o > 0 
Imho < 0 



0 

^ a 




Reh-i < 0 
Imh] > 0 

a 




0 

h-2 


lm h -2 > 0 
Reh,2 > C 





a 

b 

L 



Figure 7. Behavior of the branches h along the upper bank of R + . 

It follows immediately from (6.5) and (6.6), by choosing the initial point of integration 
to be a in (2.6), that 

(6.34) No(x)=Ni(x), x€ Aj, and Nq(x) > N-| (x), xe(a,b). 
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Similarly, by choosing the initial point of integration in (2.6) to be b, we deduce from (6.8) 
and (6.6) that 

(6.35) N 0 (x) = N 2 (x) > N 2 (x), xeA 22 , and N 2 (x) > N 0 (x), xG(a,b). 
Hence, we get from the symmetry of the surface that 

(6.36) Ai c Toi, A 2 C To 2, (—b, —a) U (q, b) G fl 2 oi. 

Moreover we have 

(6.37) N 2 (a) >N 0 (a) =N 1 (a), N 2 (b) = N 0 (b) > N, (b). 

To conclude our qualitative analysis based on Figure 7, we continue the integration beyond 
the points 1 (to the right) and a (to the left). We get that there exists di > 1 and d a > 0 
such that 

(6.38) (l,di) en 02 i and N 0 (x) = N 2 (x) > N 2 (x), x G (d a , a]. 

Indeed, if we start integrating from 1 in the positive direction, then N 0 increases and N 2 
decreases. Thus, (6.38) follows from (6.35) and (6.36) by continuity. We summarize the 
order of the branches of N along !R + in Figure 8. 


N 2 >No = T\1i Nq = 1 SJ2 >N) IN! q > N 2 > IN i No > max{lNl 2 , INI -]} 

• — • - • - • - • - • - • - 

0 x a d a a b 1 d 2 oo 

Figure 8. The result of the qualitative analysis of the branches of N along IR + . 


We note that N 2 decreases if its argument moves from d a to the left while No and N 2 
increase, see Figure 7. Thus, it is possible that there exists x a G ( 0 , d a ) such that 

(6.39) N 0 ( x a ) = N! (x a ) = N 2 (x a ). 


Moreover, if such a point exists, it is unique. In order to prove the last claim, set 

1 fa(x) dx 

dAi x = —, xG(-a,a , 

V a 2 — x 2 2 tt 

where the function f a is positive and was defined in (6.8), 

ft M 


dA 2 (x) = < 


x — b dx 


1 — X 271 


, x G (b, 1 ), 


-f-i (xl 


—x — b dx 

—, x G (— 1 ,—b). 


1 + X 271 

where the function f 2 is positive and was defined in (6.5) while and f_ 2 is negative and 
was defined in (6.10), and finally 

dx 

dA 0 (x) = ±(h 1+ —h]_)(x) ±x G [ic, ioo). 

Z 7 T 1 


The measures A-|, A 2 , and Aq are positive (for A 22 this claim follows from the discussion 
after (6.24) where one needs to recall that we deformed Aq to be as in (6.33)). It further 
follows from (6.5), (6.8), (6.10), and Privalov's lemma [31, Sec. III.2] that 


hi(z) = 


dA-i(x) 


x — z 


-Hf 


dA 0 (x) 


x — z 


and ho(z) = — 


d (A1 +A 2 )(x) 


x — z 


Since tq^z) = 23 z N 2 c (z), we deduce that the branches i G { 1 , 2 }, have the following 
global representation 

N i (z)=V i (z)-(-l) i V 0 (z)-c i , 

where V| c (z) = — Jlog|z — tldAicft), k G { 0 , 1 , 2 }. Notice that the symmetry of Of implies 
that A 2 is an even measure. Therefore, in a complex neighborhood of [—b, b] (which is the 
gap between two connected components of supp(A 2 )), the potential V 2 has the form of a 
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saddle such that on the imaginary axis iiR it is an even, concave and decreasing function. 
Hence, 


9 2 V 2 (x) 

3x 2 


> 0, x £ [—b, b], and 


The same is true for Vq, that is. 


d 2 V2(hj) 

9 y 2 


y € R. 


9 2 VqM 

3x 2 


< 0 , 


x £ R, 


and 


9 2 Vo(hj) 

3 y 2 


y £ [—c, c]. 


Thus, N2 is a convex function on [—b, b]. As 2 N1 = —N2 on [—a, a], N-| is a concave 
function on [—a, a]. Thus, inside [—a, a], the inequality Ni > N2 can be true on the 
connected set only, i.e., on the interval [—x a ,x a ], so (6.39) and the uniqueness claim are 
proved and we obtain 


(6.40) 


N 0 = Ni > N 2 on (— x a , x a ). 


Now, consider the trajectory E12 (be., the set where N -] = N2) emanating from the point 
ic into {Rez > 0 }. It cannot cross the set ilRU [ 1 ,oo] U [—00, — 1 ] as it would contradict the 
maximum principle for harmonic functions. Hence, the only possibility for to cross 
[ 0 , 1 ] is to cross it at x a . Analogous considerations lead to the conclusion that the trajectory 
of Po2 emanating from the point b to the lower half-plane arrives at x a and coincides with 
the continuation of the considered trajectory Pi 2- Thus, we have that the three subarcs of 
Pi 2 emanating from the point ic (resp. —ic) terminate at the points ioo (resp. —ioo) and 
±x a ; three subarcs of To2 emanating from the point b (resp. —b) terminate at the points 1 
(resp. — 1 ) and ±x a ; the trajectory Poi joins the points —a and +a. That is, we have shown 
that the set T has indeed the form as in Figure 6. 

The structure of P and the order of the branches {N ic}k. =0 on ® ( see Figure 8 and (6.40)) 
allow us to identify the decomposition of C\P into UTlpk like in Figure 9. It remains to 



Figure 9. Case I: The sets Qijk when Aq = [—ic,— ioo] U [ic, ioo]. 


deform the cut (6.33) to the one in Figure 3(a) while simultaneously interchanging indices 
1 and 2 in the subscripts of Oqk bounded by Ao and [—ic, —ioo] U [ic,ioo]. This finishes 
the proof of the Case I of Theorem 3. 

Now we prove Cases II and III of Theorem 3. Recall that we put p = a /(1 + a 2 )/3 in 
(2.1) and therefore the imaginary branch points ±ic annihilate at infinity. Similar to Case 
I, we can use (6.17), (6.19), (6.21), and (6.23), to show that 


Nt(z) = Vi(z)-Ci, V|(z) 


log |z —1| dAi(t), 


where At is a positive measure supported inside of At, i £ { 1 , 2 }. Then the choice of the 
additive constant in (2.7) and the ordering of the branches N along 1 R give us 


( 6 - 4 1 ) 


( 2 Vj +V 2 )(t)- 7 i 

( 2 V 2 +Vi)(t)-Y 2 


= 0, 

t € supp(Ai), 


> 0, 

t £ [-a, a]. 

j Yi =2ct +c 2 . 

= 0, 

t £ supp(A 2 ) 

\ Y2 = 2 c 2 +Cj. 

> 0, 

t £ [— 1 , —a] U [a, 1 ], 
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This type of vector-potential equilibrium problem is well-known in the theory of Hermite- 
Pade approximants and systems of measures {Ai, A2} which satisfy such equilibrium con¬ 
ditions are called Angelesco systems, [1, 27, 21]. Using the symmetry of the problem (6.41) 
with respect to the imaginary axis we can transform it to the following one. Set 

,2 


Vi (w) := 


1 


0 l0§ |w - x| 


dA] (-s/x) 


log 


1 


dA 1 (t) = (V 1 (z)+V 1 (-z ))/2 = V 1 (z) / 


. _ |z 2 -t 2 | 

where we used the change of variables w = z 2 and x = t 2 . Analogously, we put 


V 2 (w) := 


1 1 

In 7-1 dA 2 (%/x) = V 2 (z), w = z 2 

n 2 W-X 


Thus, (with the same constants y- x as in (6.41)) we have 

= 0 x € supp(Ai (v^x)) n [ 0 , 1 ], 


(6.42) 


( 2 V! + V 2 )(x) — Yi 


( 2 V 2 +V!)(x)-y 2 


^ 0 x G [ 0 , a 2 ], 

= 0 x G supp(A 2 ( 1 /x)) (T [ 0 , 1 ], 
^ 0 x G [a 2 , 1 ], 


Therefore, 


(6.43) Ni(z) =N t (w) :=Vi(w)-Ci, w = z 2 , i e { 1 , 2 }. 

Angelesco systems (6.42) for two touching intervals [ 0 , a 2 ] and [a 2 , 1 ] (for any a G ( 0 , 1 )) 
are well studied, see [24, 22, 2, 6]. In particular, one knows that the regions con¬ 
structed for Ni c (w) decompose C as in Figure 10. 



Figure 10. Domains for M. Case (I): a 2 < 1 / 2 ; Case II: a 2 = 1 / 2 ; Case III: 
a 2 > 1/2 


Thus, if we plot images of the sets £ 1 ^^ from Figure 10 (Cases II and III) under the trans¬ 
formation z = yhv, then (due to (6.43)) we obtain the corresponding sets from Figure 3 
(b,c,d). This finishes the proof Theorem 3. 

It remains to show that (6.26) defines single-valued harmonic function on the Riemann 
surface of h(-, p*), p* = \/{1 + ci 2 )/ 3 , which satisfies (6.27). This follows from the fact 
we can obtain from N(;p*) through the transformation z 2 —I w and the latter 

can be computed via (6.42) and (6.43). The partition of C by the domains II;j ^ for N is 
presented in Figure 10 [Case (I)]. Hence, the corresponding domains for N are distributed 
as in Figure 11, which finishes the proof of Theorem 1 and justifies the proof of Theorem 3. 
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Figure 11. Sets for (2.1) with p = \/(l + a 2 )/3 and a < 1 /\fi. 


7 NUTTALL-SZEGO FUNCTIONS 

This section is devoted to the construction of the Nuttall-Szego functions of Theorems 6 
and 7. Along the way we prove Propositions 4 and 5. For brevity, we shall use the 
following notation: 

9 9 

and fK„ /P := IK\ (J (a t U 00 . 

i=l i=l 


7.1 Abelian Differentials 


A holomorphic differential on IK is a differential of the form 

dQ(z) = f (z) dz, 

where f is a rational function on hi whose principal divisor is of the form 


(7.1) (f) = © f + 2 ^oo (0> + oc/ 1 ^ + 00- 2 *^ 


L 


({order of branching at e} — 1) e. 


eE {branch points of 


for some integral divisor ©f of order 2 g — 2 . It is known that such rational functions 
(integrands of holomorphic differentials) form a subspace of dimension g. Flence, there 
exists precisely one such integrand (up to a multiplicative constant) when g = 1, and 
there exist 2 linearly independent ones when g = 2 . It is further known that if f 1 and f2 
are distinct integrands of holomorphic differentials, then Df 1 and ©f , have no element 
in common. Moreover, any z £ IK has a unique complementary point, say z, such that 
©f. = z + z for some holomorphic integrand f z . When g = 2 , the integral divisors ©f in 

(7.1) are exactly the special divisors mentioned before Proposition 5. Clearly, if ©f, and 
©f 7 are any two special divisors, then 'Df t — © f , is principal. 

A point z is a Weierstrass point of a genus 2 surface IK if there exists a rational function 
on IK with a double pole at z and no other poles. Hence, z is a Weierstrass point if and only 
if z = z. To find Weierstrass points it is enough to find a two-sheeted Riemann surface 
conformally equivalent to IK. Indeed, Weierstrass points are mapped into Weierstrass 
points and the Weierstrass points of a two sheeted surface of genus 2 are precisely the 
branch points. It was shown in [10, Theorem 1.1] that IK is conformally equivalent to 

(7.2) Z 2 = A'(t) 2 - 4 A(t)( 3 t 2 -(l +a 2 )+ 3 p 2 ). 


where A is the same polynomial as in (2.1) and the correspondence between the surfaces 
is given bv : 

j t = z + 1 /h 

{ Z = 2 ( 3 t 2 — (1 + a 2 ) + 3 p 2 )/h + A'(t) 


5 One needs to replace H by —h in (2.1) to get the correct correspondence with [10, Theorem 1.1]. 
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for the uniformizing parameter t € C. It follows from (6.12) and (6.13) that 

, 2 


ttoo^) = 


1 +a 2 


and 


Z 2 too*' 2 -*'! 


= A' 


'(t(oo^)) . 

* 21 is not a Weierstrass point 


As Z(oo i2: j 0 , it is not a branch point of (7.2) and hence 00 
of SR, which is a conclusion we shall need later. 

As just explained, SR has g linearly independent holomorphic differentials. Denote by 

dO f (da 1 ,dQ 2 ) T , g = 2, 

\ dOT, g = 1, 

the column vector of such differentials or the unique differential normalized so 


dOj — 


« k 


where 6j ^ is the usual Kronecker symbol. The differentials dQ^ form a basis for the space 
of holomorphic differentials on SR. The vector of analytic integrals II defined before (3.1) 
is given by 


O(z) := 


da. 


e i 


where e-\ is the branch point of SR such that 7t( e t ) = 1 (we could have chosen any other 
point to be the initial point for integration). The significance of this vector lies in Abel's 
theorem: if T) 1 and £D 2 are integral divisors, the divisor T>-\ — £D 2 is principal if and only 
if the orders of T>i coincide and 


( 7 - 3 ) 


d(®i) = a(SD 2 ) (mod periods dd) , 


where O(SD) := ^HiCijZi) for D = Y.x n i z i- 

In what follows we shall also use differentials d£l z/vv that are the normalized (that is, 
§ da Zj ,w = 0, 1 C i C g) abelian differentials of the third kind with simple poles at z 
and w of respective residues +1 and — 1, and holomorphic otherwise. 


7.2 Cauchy Kernel on SR 

To solve the boundary value problem stated in Theorem 6, we need a (discontinuous) 
Cauchy kernel suited for our purposes. Define dC z to be the normalized third kind 
differential with three simple poles, say z,zi,Z 2 , assuming 7t _1 (z) = {z,Z],z 2 }, which 
have respective residues 2 , — 1 , — 1 . It is quite simple to see that for any fixed w we have 

dc z = 2dn Z/VV da Zl/W dn Z2 , w . 

Let y be a chain on SR possessing projective involution. For each t £ y which is not a 
branch point of SR, we shall denote by t* another point on y having the same canonical 
projection, i.e., 7t(t) = n{ t*). When t € y is a branch point of the surface, we simply set 
t* = t. 

Henceforth, fix y as above and let A be a Holder continuous function on y. For simplic¬ 
ity, we shall also assume that the set 7t(y) n Uf=i 71 1 a i) consists of finitely many points. 
Define 


A(z) := -r~r d> AdC z , 
6m 


z £ SR\7t MTtfyuQat) 


i=i 







7.3 LOGARITHMIC DISCONTINUITIES 


The function A is holomorphic in the domain of its definition. Moreover, when z->te 
a k , one of the points z- x , say z-\, approaches t* £ ajf, while z2 does not approach any 
point on y U a k (projective involution property). Hence, 


A+(t)-A-(t) = 1 


— o A(dn+ w -dn t - w )-— o A(dO t A w -dO t + 


= —o Ad£l k , 

• Y 

since, according to [38, Eq. (1.9)], one has 

dO+ w - d Cl~ w = - 27 tidQ k , z € a k . 
On the other hand, if z —> f, z-\ —> t £ ct k , and z 2 —> t* £ af c , then 


A + (t) — A - (t) = - 


-^° A ( da tw~ d£ 1 t,w)-^° A ( dn t*,w - d °u> 


J Y 

Furthermore, if z —» t £ y^ \ Uf=i a i/ ^ —► t* £ y^, while z 2 does not approach any 
point on y U Uf=i a i/ then 

A+(t )_ A -( t ,_“W + W), 

according to [38, Eq. (2.8)]. Finally, if z —> t, Z-\ —>• t £ y^ 1 , and z 2 —► t* £ y^, then 

A + (t) - A - (t) = ——. 

Thus, if we additionally require that A(t) = A(t*), then A is a holomorphic function in 
Dloc \ y such that 


= 0 . 


(A+-A-) (t) 


o AdO k , t£a k \(yU« 3 _ k ), 
- Y 

A(t), t £ y \ (Ji =1 a i- 


It also can readily be verified that 

(7.4) A(z) + A(zi) + A(z 2 ) = 0 on % 

since the cycles y and oq possess projective involutions. 


7.3 Logarithmic Discontinuities 


It is known that the continuity of A^, in fact. Holder continuity, depends on the Holder 
continuity of A only locally. That is, if A is Holder continuous on some open subarc of y, 
so are the traces on this subarc irrespective of the smoothness of A on the remaining 
part of y. Of particular importance for us are the logarithmic discontinuities at the branch 
points of 1H. 

More precisely, let e be a branch point of 91 and it be its circular neighborhood. Given 
a cycle on 91 that possesses a projective involution and passes through e, denote by y its 
part that belongs to it. In this case 7t(y) is a Jordan arc ending at e = 7t(e). Recall [20, Sec. 
1.8.6] that 


( 7 - 5 ) 


1 

27 ti 


alog(t—e) dt 
7t(y) G 


v+ 


t —Z 


= ± 


ae 


±V7ri 


log(z-e) 


2 i sin(v7r) 


z-e 


+ 0((z — e)- v ) 


as z —> e, where cdog(z — e), a £ R, and (z — e) v , 7/O, are holomorphic outside of 7 t(y), 
log(t — e) is any continuous determination of the logarithm on 7 t(y), (t — ej v " is the trace 
of [z— e) v on the positive side of 7 t(y), and the choice of the sign in ± is determined by 
the orientation of 7 t(y): + if 7 t(y) is oriented away from e and — if it is oriented towards e. 
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Assuming that 91 has double branching at e, define u(z) := ±(z — ej 1 / 2 , where the 
choice of the sign depends on whether z lies to the left (+) or to the right (—) of y. Set 


dC 2 (t) := 1 


1 +3 


u z 


dt 


u(t) J t — z 

The differential dC 2 is holomorphic in U\{z, z*j with residues +2 at z and —1 at z*, where 
z* is the unique point in it which is different from z and has the same natural projection. 
Clearly, the difference dC 2 — dC z is a holomorphic differential in it. Therefore, if Aft) is 
of the form ixlogft — e) + {Holder continuous part} satisfying A(t) = Aft*), then 


(7-6) 


A(z) = -n 

b7tl 

1 

2 ™ 


a log ft - e) de 2 (t) + 0 ( 1 ) 


. , , ufz) 

alog(t-e) 

tt(y) vt —e 


■^^ + 0 ( 1 ) = ±|log(z-e) + 0 ( 1 ) 


as z —> e S qA by ( 7 . 5 ) applied with v = 1 /2 (observe that the second equality above is 
valid independently of the orientation of 7t(y)). 

Assume now that 91 has triple branching at e. Since y possesses projective involution 
there is exactly one component of il\y that lies schlicht over 7 i(il\y). Orient y so that 
this component lies to the right when y is traversed in the positive direction (it is bordered 
by y _ ). Define 


de 2 (t) := 


u z 


+ 


u 2 (zl 


u(t) u 2 (t);t 


dt 


where ufz) is equal to the lift of a fixed determination of (z — e ) 1//3 to the "schlicht" 
component of il and then is continued holomorphically to the whole it. Since the values 
of u at the points with the same natural projection differ by e ±2 ™/ 3 , one has that dC 2 
is a holomorphic differential in ll\ 7t _1 (z), with 3 simple poles at the elements of 7T _1 (z) 
having residues +2 at z and —1 at the two elements of 7+ 1 (z) \{z}. Again, dC 2 — dC 2 is 
a holomorphic differential in it and therefore 

alog(t— e) dC z (t) + 0 ( 1 ) 

‘ T 


A(z) = 


67 ti 

\/3i 

6ni 


tt(y) 


alogft — e) 


1 ±7ti/3 u ( z ) 1 ±27ti/3_ - , _ 

(t-e) 1 /3+^ e (t — e) 2 /3+y t — z' 


u 2 (z 


(z) A dt 


where the sign + is used when 7t(y) is oriented towards e and the sign — is used if it is 
oriented away from e, and we used the fact that 


dt 

dt 

dt +' / 3ir ±7Ti/3 

dt 

.yU(t) . 

7 i(y) (t — e ) 1 / 3 

(t-e) 1 /3+ J 

7t(y) (t — e) 1 / 3 


with the same choice of signs. Then ( 7 . 5 ) applied with v = 1/3 and v = 2/3 gives that 


A(z) = 




--log(z-e 


(7-7) 


oc, 


0 ( 1 ) + ylog(z —e 



+ 0 ( 1 ) 


where the determination of (z — e ) 1 / 3 is precisely the one of ufz) within the component 
of it that lies to the right of y. 


7.4 Szego Function S p 

Let A be the chain introduced in Theorem 6. We assume that the cycles forming A and 
separating 93^°' and 9V i: , i £ [1,2}, are oriented so that 93 ' 0 lies on the left (positive) side 
as each of the cycles is traversed in the positive direction, while the cycle separating 93 ' 11 
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and 91 : 2 is oriented so that 91 : 1 1 remains on the left when it is traversed in the positive 
direction. 

To define a density whose Cauchy transform on A we shall need, let us introduce 
auxiliary "trigonometric" weights. Namely, 


(7-8) 


W] (z) := Vz 2 — a 2 , 

w 2 (z) := \J[z 2 - 1 )(z 2 -b 2 )(z 2 + c 2 ), 


in Case I with branch cuts along A], A2 U Ao, respectively (all are positive for positive 
reals large enough), and 


(7-9) 


W] (z) := Vz 2 — b 2 , 

W 2 (z) := ^/{z 2 - l)(z 2 - a 2 ). 


in Cases II and III with branch cuts along Ai and A 2 , respectively (all are positive for 
positive reals large enough). We set 

r 


(7.10) 


A p (t) := < 


log (piwf)(t). 

t e A° v 

log(-p 2 w+)(t). 

t e A° 1V 

log (p 2 W+)(t), 

t G A° 2 , 

log (p^J/ptW!)^), 

t G A° 0 , 


t € A, 


where we choose continuous branches of the logarithms, which is possible as the functions 
Pi and the trigonometric weights are non-vanishing except at the endpoints. Put 


(7.11) 


S p (z) :=exp{A p (z)}, c p := -^7 o A p dQ. 


Then S p is a holomorphic and non-vanishing function on 91 a \ A with continuous traces 
except at the branch points that satisfy 

f exp { 27 ti(c p ) k } on a k 
exp{A p j on A. 

Moreover, S p satisfies (2.10) with ® replaced by S p as follows from (7.4), and, excluding 


(7.12) 


S+ - c- 


(7-13) 


/yfl, one has 



' jSp°^(z) ~ z — e | -1 / 4 

as 

z —t e G {±1, ±q, ±b}. 

]Sp 1 ] (z)| ~ z —el 1 / 4 

as 

z — y e G {ic]}, 

|Sp 1 ] (z) ~ z — e | -1 / 4 

as 

z — > e G {±ic}. 

. ]Sp 2 ) (z) ~ z —e | 1//4 

as 

z-teG (±1, ±e 2 ,±ic}. 


by (7.6), where e-\ = a and e 2 = b in Case I while ei = b and e 2 = a in Cases II and 
III, and the points ±ic appear only in Case I. When e-\ = e 2 = 1 /V' 2 , we need to use 
(7.7) instead of (7.6). Now there are two cycles that pass through a branch point: one 
that separates and 91 11 ' ( 91 * 1 ’ lies to the right of this cycle), and one that separates 
91 *°* and 91 2 ( 9 T :2 lies to the right of this cycle). Combing the contributions from the 
integrals over both cycles we get that 


Sp°^(z)| - |z — ej as z->e€{±l /y/l}, 

Sp l) (z)| ~ |z— el 1 / 6 as z->ee{±l 1 € { 1 , 2 }. 


(7-i4) 
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7.5 Riemann Theta Function 


The theta function associated with B is an entire transcendental function of g complex 
variables defined by 

0 (u) := ^ exp < 7rih T Bh + 27 iin T u >, u G C 9 . 

ftezs 

As shown by Riemann, the symmetry of B and positive definiteness of its imaginary part 
ensures the convergence of the series for any u. It is known that 

d (w) + K, g = 2 , 

K, g = 1 , 


(7-i5) 


0 (u) = 0 <=> u = 


I mod periods dd 


for some w € 91, where the equivalence = was defined in (3.3) and K is the vector of 
Riemann constants defined by 

f (K) k := (Mick - 1) /2 -#„ 3 _ k ® k dD 3 _ k , g = 2, 

1 K := (B—1)/2, g = l. 

Let 2 ) be an integral divisor of order g. Then 0 (d(z) — dfD) — K) is a multi-valued 
holomorphic function on if 2 ) is not special and is identically zero otherwise. Indeed, 
in the later case 0 ( 2 )) = 0 (z + z) + j z + Brh z , j z , rft z G Z 9 and therefore 


0 


(d( 


z) — 0(2)) — K) = 


'j =0 d(£)— e -™™I B ™z- 27ti ™-I(^(z)+k) _q 


by (7.15), where we used the fact 0 (—u) = 0 (u) and the periodicity property of theta 
functions: 


(7.16) 


0 ^u +j + Brri^ = exp < — 7 rim T Brh —27tim T u >0(u), j,m 6 Z 9 . 


Set v k := d(oo( k ’ — oo* 0 ^), k G { 0 , 1 , 2 }, and define 


(7-17) 


®n,k(z) := 


0 ^O(z) -O(goo (2) ) - K-v k -c p -x n -Bq T 
0 (d(z)-d(goo( 2 )) -K) 


(recall that 200 2 is not a special divisor, see ( 7 . 2 ) and after, and therefore 0 n<k is always 
well defined), where x n , y n G [0,1) 9 are such that 


(7.18) 


Xn + Byn = n (u3 + Bt) (^mod periods dO^j . 


If the numerator is not identically zero, 0 ll (c is a multiplicatively multi-valued mero- 
morphic function on 9H with g simple poles at 00 :21 , g zeros that we shall describe by an 
integral divisor 2 ) n ,k, an d otherwise non-vanishing and finite (it is also possible that zeros 
could cancel poles). Moreover, it follows from (3.1) and (7.16) that 0 n<k is meromorphic 
and single-valued in thd and 

(7.19) 0+ k = 0- k exp{-27ti(v k + Cp+Xn + Byn) i } 


on a;. 


Together with the functions 0 n ,k we shall need two more auxiliary theta functions. We 


set 


0 k (z) = 


e(n(z)-n( OO + (g — 1 )w) — 

0 (a(z) - d(oo( k ) + (g - 1 )w) - K) ' 


kG {0,1,2}, 


for any fixed w such that w 7^ c»f k ' for all k G { 0 , 1 , 2 }. Clearly, 0 o = 1 and 0 i, i G { 1 , 2 }, 
is a meromorphic function in D\oc, with a simple zero at 00' 0: , a simple pole at 00' 1 , and 
otherwise non-vanishing and finite. Moreover, its traces on a k satisfy 

(7.20) 0 += 0 Texp{ 27 ti(v i ) k } on a k . 
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7.6 Jacobi Inversion Problem 

It follows from (7.15) and (7.18) that the divisor D n ^ is a solution of the Jacobi inversion 
problem 

(7.21) 0 ( 1 )^) = D(goo (21 )+Vk+ c p +n(d 3 + Bf) ^mod periods dO^ . 

Clearly, '11,1,0 are precisely the divisors 'D n defined in (3.2). 

As we already mentioned after (3.2), any Jacobi inversion problem is uniquely solvable 
on surfaces of genus 1. Thus, are well defined in this case. Since 

(7.22) 0(11^1;) = + (!) + Bt (mod periods dO^ , 

Proposition 2 and the unique solvability property imply that 

®2m,k = ®0,k 7^ ®1,k = ®2m+l,k 

for all m ^ 0 . This finishes the proof of Proposition 4. Observe also that 

(7.23) 0(2)^) = Cl(T) n j — 00^ 1 + oo^ k) ) ((mod periods dQ^ . 

The latter equivalence together with the unique solvability property imply that 

T> n/ j=oo^* for some jG { 0 , 1 , 2 } => ^D n ,k = oo^ k * for all kG { 0 , 1 , 2 }. 

Hence, if T> n ,0 7^ oo^ 0 *, then T>n.,x 7^ oo* 1 ', i G { 1 , 2 }, which will be important in what 
follows. 

As to Proposition 5, observe that its first claim follows from (7.22) and Proposition 2 
as fi assumes equivalent values at special divisors by Abel's theorem (7.3). To prove the 
second claim, assume to the contrary that all the limit points of {D tli o } ne lrx - are either 
special or of the form 00 ^ +w for some w € 91 . Using compactness, we always can 
select a subsequence N1 such that 

t 2P as IN - ] 3 n —k 00 

for some divisors TL simultaneously for all j G { 0 , 1 , 2 , 3 }. According to our assumption, 
one has 

(7.24) £>(2P) = D(oo (0) +Wj), 

where Wj = 60 ; 01 when ‘D ] is special. Observe also that continuity of Q and (7.22) imply 
that 

(7.25) Q(D^ + 1 ) = O(O^) + u 3 + Bt (mod periods dfl() 

for ) G { 0 , 1 , 2 }. By combining (7.25) for different j and then recalling (7.24), we deduce that 

a(2wj) = n(wj + i +wj_i) ((mod periods d£l(), 

j G { 1 , 2 }. The latter equivalence and Abel's theorem imply that 2 wq and W3 +W] are 
special divisors, and therefore necessarily W3 = W]. The above conclusion and (7.25) 
imply that 

2 cu + 2 Bt = 0 ((mod periods d£l(j , 

which is impossible since 2 u 3 = ( 2 cu, 4(1 — tu)) and 2 cu G ( 1 , 2 ) by Proposition 2. This 
contradiction proves the second claim of Proposition 5. 

Finally, observe that the subsequence IN, is such that no limit point of { is 

special or of the form 00 ■ 1 + w for i G { 1 , 2 }. Indeed, otherwise there would exist a limit 
point D of one those sequences such that 

O(D) =0(oo (i) +w). 
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where w = ob^) if D is special. Choose a limit point D* of {D n ,o} along the same 
subsequence. The continuity of Cl and (7.23) then would imply that 

0(2)*) = 0(D — 00^ +oo (0) ) = d(oo (0) + w) (^mod periods dD^j, 

which means that either CD* is special or contains 00 As both options are impossible 
by the choice of N*, such a limit point D does not exist. 

7.7 Nuttall-Szego Functions 

Finally, we are ready to prove Theorem 6. Let y £ ]R tJ . Then it follows from (3.1) that 
Sg(z) :=exp{- 27 tiy T d(z)}, z e 9 t a/P , 
is a holomorphic function on D\ a p with continuous traces that satisfy 

, s + = s-J exp ^ 2rti ( B ^^ on 

g g \ exp {— 27 ri(y) i } on ( 3 ^ 

Moreover, (7.4) implies that (2.10) holds with CD replaced by Sg. Finally, set 

W(z^ 0 ') = 1 and W(z^ x) ) = w;(z), ie{l, 2 }, 

where the functions W] and W2 are defined by (7.8) or (7.9), depending on the considered 
case. The Nuttall-Szeg6 functions we are looking for are defined by 

(7-27) Vn,k := O^SpSg^n^OkW" 1 , 

where y n was defined in (7.18) and we assume n is such that the corresponding Jacobi 
inversion problem (7.21) is uniquely solvable. Indeed, it can readily be checked, using 
(2.11), (7.12), (7.18), (7.19), (7.20), and (7.26), that each is holomorphic in 9 T \ A and 
its traces satisfy 

(WT n , k ) + = (WT n , k )-exp{A p } on A. 

Using the definition of W and (7.10), we see that the pullbacks of k solve the boundary 
value problem (3.6). Further, the functions T n k have the local behavior around the branch 
points of as stated in (3.5) by the very definition of W and (7.13) or (7.14). Finally, the 
zero/pole multi-set of U TX<k in 91 \ A is described by the divisor 

(7.28) ©nA + (n+ 1) ^00^ + oo^ 2 '^ — 00^ — (2n — 1 )oo^ 0 ' 

as follows from the properties of <D n , 0 n k , 0 k , and W. In particular, the functions 
satisfying conditions of Theorem 6 are given by T n := U tli o- 

It only remains to show the uniqueness of dV /k . To this end, consider dV k/T, where 
¥ satisfies (3.5), (3.6), and (7.28) with B U(k replaced by any other integral divisor, say T>. 
Then the analytic continuation property implies that T' n k/T' is a rational function on 9 L 
whose divisor is given by ® n(k — T>. As © n/k is not special, T> = D n x, that is H'^k/'T i s a 
constant. 

7.8 Auxiliary Observation 

In Section 9, we shall need the following observation. Assuming we are in Case II, let 
®n,k(p) be the solution of (7.21) for the weights pi and p2 and 2)^0/p) be the solution 
of (7.21) for the weights 1 /pi and 1/P2- Then it can easily be seen from (7.21), (7.10), and 
(2.15) that 


n(2) n ,k(p) + 2Vk0/p)) =2a(oo( 2 >)+2v k + 2c 1 = 2n(T> rL/k (l)). 


7.9 ASYMPTOTIC BEHAVIOR 


Moreover, it also follows from (2.15) and (3.6) that (S-| S LJ ' ii 0 tl ^ 0 [ c )“ is a rational function 
on 93 with the divisor equal to the sum of the divisors (7.1) and 2D rl/ i c ( 1 ) — 2oo k 1 . As the 
divisor (7.1) is a divisor of a rational function, we get from Abel's theorem (7.3) that 

H(2)n,k(p) + »n,k(1/p)) =2D(0O^). 

In particular, it follows from the unique solvability of the Jacobi inversion problem that 
(7.29) ®n,k(p) =OO tk) I>n,k(Vp) =OO tk) . 

That is, n g N*(p) if and only if n g N*(1/p). 


7.9 Asymptotic Behavior 


Below we prove Theorem 7. In fact, we shall show that a more general statement is true. 
Given IK, and e > 0 small, there exists a constant 1 < C £ (IK*) < 00 such that 


( 7 - 30 ) 


|M'n,kl<C e (N*)|®r in 93 £ , 
l^kl^ CefN*)- 1 \(D\ n in 4 k) , 


for all n g N* large, where il[ k * := 93 * k ) 0K 1 ({|z| > 1 /e}). 

It follows from the continuity of the boundary values of S p as well as from (7.13) that 
|S p | is bounded from above and away from zero in 93 £ (that is, including the boundary 
values on A and the a-cycles). Moreover, since the image of 93 under Cl is bounded 
and y n g [ 0 , 1 ) 9 by the very definition, see (7.18), the moduli Sg n | are uniformly bounded 
from above and away from zero. Thus, we only need to estimate the function ©n,k©kW _1 . 
It follows from (7.28) that the zero/pole multi-set of 0 n,k©kW -1 in 93 « \ A is equal to 

k + 00 ■ 0 * + 00 1 + (x/ 2 * — 00' k 1 and therefore these functions are holomorphic there 

(k) 

as well as non-vanishing in il £ for all £ > 0 small and n large by the very choice of IK,. 
Hence, 


(7-3i) 


©n,k©kW C Cn. /£ ( 1 N*) in 93 £ , 
©^©kW-’l ^ C n , £ (NJ- 1 in 4 k) - 


Thus, if we show that { 0 n< k 0 kW _1 } is a normal family and every limit point satisfies 
the estimates as in (7.31), a standard compactness argument will finish the proof of (7.30). 
The only varying part of this family is given by 

(7.32) 0 (fi(z) -Q(goo (2) ) -K-Vk-Cp-Xn-Bljn) . 

Boundedness of its image follows from the continuity of 0 ( ), boundedness of the image 
of d(z) for z g 93 a< p, and the fact that x n ,y n g [ 0,1 ) 9 . Thus, { 0 tv ^ 0 ^W _ 1 } is a normal 
family and its elements are indexed by pairs (xnv'yn) £ [o, I) 29 . The continuity of 0 (-) 
and (7.17) imply that any limit point of the functions in (7.32) has the same form with 
(Mutjn) replaced by a limit point (x, yj g [ 0,1 ] 29 . Denote by D the integral divisor of 
order g describing the zeros of this limit point. Then 

0 ( 2 )) = 0 (goc/ 2) ) + Vk + c p +x + By (^mod periods dtl^ . 


The continuity of li implies that the right-hand side of the equivalence above is a limit 
point of the right-hand sides in (7.21). By the definition of IK,, it determines T) uniquely 
as it is a limit point of {©n,k}n.eN* • Moreover, all such limit points are uniformly bounded 
away from containing 00and therefore from containing an element from U £ k for all 
£ > 0 small enough, which is sufficient to get the required lower estimate. 
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8 RIEMANN-HILBERT ANALYSIS: CASE I 


8.i Global Lenses 



( 1 

0 

° ^ 

/ 1 

0 

0 \ 

(8.1) 

G1 (u) := 0 

1 

0 and G 2 (v) := 

0 

1 

v 


V 0 

u 

1/ 

1 0 

0 

1/ 


Observe that Gt(u) 1 = Gt(—u), i € { 1 , 2 }. Moreover, with J as in (5.5), 

G^ 1 (u)J(x,y)Gi (u) = J(x + uy,y) 

G? 1 (v)J(x,y)G 2 (v) = J(x,y + vx). 



That is, conjugating J by G^ allows one to modify the i-th variable of the matrix J using 
the other variable. 

Let, as before, ±ic, c > 0 , be the projections of the branch points of 91 and V be the cut 
along which the sheets 91 * 1 * and 91 : 21 are glued to each other. Pick arcs Aqi ± and Aq 2/ all 
oriented from —ic to ic, as in Figure 12. These arcs delimit three bounded domains that 
we label from left to right as Oq, 01, and 0 2 . 



Figure 12 . The lens Is and the domains Oq, Oi, and Oi. 


Set S := Y outside of the closure of Oo U Oi U 0 2 and inside put 

f Gt (—pi /p 2 ) in Oo, 

(8.3) S:=Y< G 2 (—p 2 /pi }Gi (pi/p 2 ) in Oi, 

G 2 (—P2/P1) in 0 2 , 

where the pi are the weight functions in (4.1). This way S has jumps on A01 ±, Ao 2 , and 
Aq . In particular, the jump on Aq is equal to 


Jo := Gi (—Pi /P 2 )G 2 (p 2 /pi )Gi (—Pi /p 2 ) 


10 0 \ 

0 0 p 2 /pi 

0 pi/p 2 0 J 


For brevity, let us write 


(8.4) 

Put I s 


p 2 [zj := j 

[— 1 , 1 ] U Aqi+ U Aqi_ U Aq 2 


p 2 (z], Re(z] > 0 , 

-p 2 (z), Re(z] < 0 . 

U Aq. Then S solves RHP-S: 


(a) S is analytic in € \ Is and lim z _>oo S(z)diag (z 2 n ,z n ,z TL ) = I; 
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(b) S has continuous traces on 11 := Is \{±1, ±a, ±ic} that satisfy 


' J(0,p$) 

J(P1,0) 

s+ = s_ < Jo, 

Gi (P 1 /P 2 ) 

G 2 (-p2/Pl. 

(c) S satisfies RHP-Y(c) (see Section 5 ) with [—1,1] replaced by Ig. 

If RHP-S is solvable, then so is RHP-Y, and the solutions are connected via ( 8 . 3 ). 


on (-1,-a) U (a, 1), 
on (-a, a) \ { 0 }, 
on A 0 , 
on A 0 i±, 
on A 02 ; 


8.2 Local Lenses 

For the next step we introduce additional arcs A]± and systems of two arcs A 2 ± as on 
Figure 13 , all oriented from left to right. Denote further by O i ± the domains bounded by 



Figure 13. The lens Ispecifically, the arcs A1-1- and systems of arcs A 2 ±. The 
black curves constitute the system In . 


(—a, a) and the arcs Ai±, respectively, and by 0 2 ± the open sets bounded by A 2 ± and 
(—1,— b) U (b, 1). Set 

(8.5) Z:=SLf in O i± , 

i G { 1 , 2 }, where 



/ 1 

0 

0 1 

1 / 1 

0 

0 \ 

( 8 . 6 ) 

Li := 1 /Pi 

1 

0 1 

and L 2 := 0 

1 

0 


V 0 

0 

1 ) 

V 1/P2 

0 

W 


It can readily be checked that Z + = Z_Jt on At, 1 G {1,2}, see ( 2 . 8 ), where 

(8.7) Ji := 


( 0 

pi 

0 1 

1 1 

( 0 

0 

p*2 ) 

i / pi 

0 

0 1 

and J 2 := 

0 

1 

0 

V 0 

0 

') 

1 1 

V -1/P2 

0 

0 J 


Put Iz := Is U Ai + U Ai _ U A 2+ U A 2 _. Then Z solves RHP-Z: 

(a) Z is analytic in C \ Iz and lim^oo Z(z)diag (z 2 n ,z _n , z _n ) = I; 

(b) Z has continuous traces on each side of Ih := Iz \ {± 1, ±a, ±b, ±ic} that satisfy 


Z+ 



J k on A°\{0}, ke {0,1,2}, 
Jz on ^z \ (Ao U Ai U A 2 ), 


where 


J( 0 , p* 2 ) on (—b, —a) U (a,b), 

Gi(pi/p 2 ) on Aoi±, 

G 2 (—p 2 /pi) on A 02 , 

Li on A i± , i e {1,2}; 


( 8 . 8 ) 


Jz := 
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(c) Z satisfies RHP-S(c) with Is replaced by Iz- 

As before, if RHP-Z is solvable, then so is RHP-S, and the solutions are connected via 
(8.5). 

8.3 Global Parametrix 


Let In := Ao U Ai U A2, see Figure 13. In this section we are seeking the solution of the 
following Riemann-Hilbert problem (RHP-N): 

(a) N is analytic in C \ In and lim z _^ 0O N (z)diag (z _2n , z n , z n ) = I; 

(b) N has continuous traces on each side of In \{±l,±a, ±b,±ic} that satisfy N + = 
N_J k on A°,ke { 0 , 1 , 2 }. 


We solve RHP-N only along the subsequence N* defined in Proposition 5. Let 'Fn k be 
the n-th Nuttall-Szego functions constructed in (7.27). As the functions 'P T1 k satisfy (3.6) 
with the zero/pole sets described by (7.28), it readily follows that a solution of RHP-N is 
given by 


(8.9) 


N := C r 


V. 

V. 

V V. 


(0) 

n,0 

(0) 

n,1 

(0) 

n,2 


/Cl) 

n,0 

/I) 

‘ n,1 
T n,2 


K 


2) 

n,0 

2 ) 

n,1 

2 ) 

n,2 


=: C _TL MD n , 


where n £ N„ C n is a diagonal matrix of constants chosen to fulfill the normalization 
condition in RHP-N (a), C := diag(Co, Ci, C2), see (2.13), and D := diag 1 \ ® 12 ^. 

Observe that det(N) is an entire function in C \{± 1 , ±a, ±b, ±ic} since the determinants 
of the jump matrices in RHP-N(b) are all equal to 1. Moreover, the normalization at 
infinity implies that det(N)(oo) = 1 and therefore det(N) is a rational function. It also 
follows from (3.5) that near any of the points in {±1, ±a, ±b, ±ic) the entries of one of the 
columns of N are bounded and the other entries behave like 0 (|z — e| _1 / ' 4 ). Hence, 


det(N)(z) = 0 (jz — e| 1//2 ^) as z —t {± 1 ,±a,±b,±ic}. 


and therefore det(N) = 1 . 


8.4 Local Parametrices 


Denote by U e , e £ {±l,±a, ±b,±ic}, a disk centered at e of small enough radius so that 
I e := U e n Iz consists of disjoint, except at e, analytic arcs. We are seeking a solution of 
the following RHP-P e : 

(a) P e is analytic in U e \ I e ; 

(b) P e has continuous traces on each side of 1 ° that satisfy RHP-Z(b) within U e ; 

(c) P e is either bounded or has the behavior near e within U e described by RHP-Z(c); 

(d) P e = M (I + 0(1 /n)) D n uniformly on 0 U e \ Iz/ where M and D are defined by 
(8.9). 


We solve each RHP-P e only for n e IN*. For these indices the above problem is well- 
posed as det(N) = 1 and therefore N -1 is an analytic matrix function in C \ In- In fact, 
the solution does not depend on the actual choice of IN*, however, the term 0 ( 1 /n) in 
RHP-P e (d) may depend on the choice of this subsequence. 

In solving RHP-P e , it will be convenient to use the notation 03 = and intro¬ 

duce the following matrix transformations T] A, T2A, and T3A given by 


/I 

0 

° 1 


/ 

Al ii 

0 

[All2 N 


f All 1 

Ah 2 

° ^ 

0 

[Ain 

[A] 12 

/ 


0 

1 

0 

, and 

Al 21 

Al 22 

0 

V° 

[A] 21 

Al 22 ) 


l 

Al 21 

0 

Al 2 2 / 


l 0 

0 

1 1 
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respectively, where [A]^ is the (i,k)-th entry of A. Observe that Tj(AB) = T j (A) T j (B) 
for all 2 x 2 matrices A, B. 


8.4.1 RHP-P ±1 

In [25], a 2 x 2 matrix function was constructed out of Bessel and Hankel functions that 
solves RHP-V: 

(a) 4 f is holomorphic in € \ ( 1 + U I_ U (—00, 0 ]), where I± := {C : arg(C) = ± 27 t/ 3 }, all 
the rays are oriented towards the origin; 

(b) 4 ' has continuous traces on I + U I_ U (—00, 0 ) that satisfy 

f / 1 Ch 


4 ?+ = W- 


1 1 

0 1 

-1 0 


on I ± , 
on (—oo,0); 


+ 0(r 1/2 ))exp{2C 1/2 (J3} 


(c) ¥(C) = 0 (log|CI)asC^O; 

(d) 4 ' has the following behavior near 00 

■to= p^p 72 a (; ;; 

uniformly in C \ ( 1 + U I_ U (—00, 0 ]). 

Furthermore, 03 4 '03 solves the same R-H problem only with the reversed orientation of 
all the rays. Notice also that RHP- 4 '(d) should be replaced in this case by 

r ff 3 / 2 ^/l -P 

x/ll-i 1 

To carry 4 / from the C-plane to U e , e £ {± 1 }, we need to introduce local conformal maps. 
To this end, set 

, 1 <•* 


0 3 4 / (C )03 = (2tt£ 1/2 ) 


(l + o(r 1/2 ))exp{ 2 C 1 / 2 0 3 }. 


9 e (Z) := 


(h 0 -h 2 ) (t)dt = jlog (W 0) (z)/® (2) (z)) , 


z € lf e \ [— 1 , 1 ], where the second equality follows from (2.5) and (2.9). As mentioned 
after (6.4), g e has purely imaginary traces on (— 1 , 1 ) fl U e that differ by a sign. Moreover, 
since g e vanishes at e as a square root, g 2 is conformal in li e . Furthermore, 

g 2 ({x: sgn(e)x > 1 }(T U e ) C {z : z > 0 }, 
gi((—i, i) n u e ) c {z: z < 0}, 
g 2 (A 2 ±nUe) c {z : sgn(e)arg(z) = ± 2 tt/ 3 }. 

Indeed, the first property follows from (2.4) and (6.11) while the second is a consequence of 
the fact that g e has purely imaginary traces there. The last property is the requirement we 

1 j2 

impose on the arcs A 2 ±. Choosing the branch of g e ^ which is positive on {x : sgn(e)x > 
l}n U e , one has on (— 1 , 1) (T U e that 


1/2 , ,. 1/2 
9 e+ =sgn(e)ig e l. 


(8.10) 

Now, it can readily be verified that the matrix function 

(8.11) Pe := E e T 2 V e (n 2 g 2 / 4 ) W e , 

satisfies RHP-P e (a,b,c) for any holomorphic matrix function E e , where 4 '] 4 ', 4 , _i := 

0 3 4 / ct 3 , and 

(8.12) W e \= A\ag[(® [0] /^ 2 , ( 0 (1) ) n , (O (0) O (2) ) n/2 y/p^ . 

Moreover, one has on 3 U e that 

/-o 3 /2 fm ,„ r <r 3 /2lf 1 sgn(e)P| ^ (i + o(l/n))D n . 


(8.13) Pe=E e T 2 P2 


(Ttnge 


\fl \ s g n ( e )i 


1 
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Thus, it only remains to choose E e so that RHP-P e (d) is fulfilled. Direct computation 
using RHP-N(b) and (8.10) shows that the matrix 

(8,4) E, := MT 2 ( P y> /2 (m9 «)«/ 2 _L s grt( e J')) 

is holomorphic in U e \{e}. As the entries of M behave like |z — e|~ 1 0 as z — > e by (3.5), 
the entries of E e can have at most a square root singularity there. Thus, E e is holomorphic 
throughout U e as desired. 


8.4.2 RHP-P ±a 

To solve RHP-P-ta, we again use the matrix ML In this case we have an additional 
complication coming from the jumps on (—b, — a) U (a, b). To circumvent it, we shall need 
the following fact about the matrix Mb 

(8.15) m 11 ( 0 =I 0 ( 2 C 1/2 ) and [M^] 21 (C) = 2 7 riC 1/2 I , 0 C 2 C 1/2 ) 

within |arg(z)| < 27 t/ 3 , where Io is the modified Bessel function of order 0 . Observe that 
both functions above are in fact entire in the whole complex plane. 

Define 


(8.16) 


9e(z) := j 


(h 0 - h,) (t)dt = - log (O' 0 ' (z)/O n 1 (z)) , 


z £ U e \ [—a, a]. As mentioned after (6.7), g e has purely imaginary traces on (—a, a) n U e 
that differ by a sign. Moreover, since g e vanishes at e as a square root, g 2 is conformal in 
U e . Furthermore, 


g 2 ({x: sgn(e)x > a}nU e ) 

C 

{z : z > 0}, 

g 2 ((-a, a)nU e ) 

C 

{z : z < 0}, 

g 2 (A 1± nU e ) 

C 

{z: sgn(e)arg(z) = ± 27 t/ 3 }. 


where the first property follows from (6.6) and (6.9), the second is a consequence of the 
fact that g e has purely imaginary traces there, and the third is a requirement we impose 
on the arcs A] ±. Choosing the branch of g e ' which is positive on {x : sgn(e)x > a}IT U e , 
we see that (8.10) holds on (—a, a) IT U e . 

We further define 

Fni(z) := P2 ( z ) (O ( - 2 - ) (z)) _3rL/2 I 0 (ng e (z)) 23 H log f=a, 

Fn2(z) := P2(z)(0 (2) (z)) _Dn/2 7tinI^(ng e (z)) 2 i H log|5^, 

where Io is the modified Bessel function of order o, see (8.15). The above functions are 
holomorphic in U a \ [a, b] and 

F ni+ -Fni- = mil (n 2 g 2 / 4 )p 2 (o' 2 ))“ 3n/2 
on (a,b). According to [16, Eq. 10.40.5], 

(ug a ) 1/2 F n i = (e n 9 “ 0 (l) + e- n 9 - 0 (l)) (0' 2 ')“ 3n/2 

/®c°A n /q,nr " 


0 ( 2 ) 


0 ( 1 ) + 


0 ( 2 ) 


0 ( 1 ) = 0(1 


uniformly on 3 U a by (8.16) and since 3 U a C <+201 by Theorem 3 and the choice of the 
radius of U Q . Clearly, o( 1 ) in the above equality is geometric. Moreover, a completely 
analogous estimate holds for F n2 . Given these two functions, we can set 

/ ffln (n 2 g 2 / 4 ) [M / ] 12 (n 2 g 2 / 4 ) F n1 \ 

V a := m 21 (n 2 g 2 / 4 ) m 22 (n 2 g 2 a / 4 ) F n2 . 

V 0 0 1 / 
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Then the solution of RHP-P a is given by P a := E a 4 ' a W a , where 

| W a := diag((o( 0 >(I)( 1 l) n/2 /VPT,( ( I )(0)( l )(1) ) n/ VpT,(® (2) ) n ) 

E q := MT 3 (p^ / 2 {nng a )^ /2 ^^. . 

To verify RHP-P a (d), observe that 

T 3 (7tng a ) CJ3/2 jY a T3(c> (0 V® (1) )~ 1X<T3/2 = I + O(Vn) 

by the asymptotic properties of F n i. The matrix P_ a solving RHP-P_ a can be constructed 
analogously. 


8.4.3 RHP-P ±b 

In [13], a 2 x 2 matrix was constructed out of Airy functions that solves RHP-<L>: 

(a) <t> is holomorphic in C \ ( 1 + U I_ U (—00,00)), where the real line is oriented from 
left to right; 

(b) ® has continuous traces on I + U I_ U (—00, 0 ) U ( 0 ,00) that satisfy 


® + = ®_ < 


1 0 
1 1 

0 1 

-1 0 

1 1 
0 1 

(c) <t> is bounded around the origin; 

(d) <t> has the following behavior near 00: 


on I±, 
on (—oo,0), 
on (0,00); 


®(C) = C a3/4 -b 


1 n 


I + O 


(v«)) 


exp 




\/2 \i 1, 

uniformly in C \ ( 1 + U I_ U (—00,00)). 

Again, cr 3 ® cr 3 solves the same R-H problem only with the reversed orientation of all the 
rays. As in the case of Mt, 

f r3/2 


<r 3 ®(Off3 = r« /4 T (_t ,‘) (l + O (r 3/2 )) exp |-?c 

To map ® into U e , e C {±b}, define 

(h 0 -h 2 ) (t)dt = ^ log (® (0) (z)/® (2) (z)) , 


^3 


1 rZ 

g e (z) := -- 


z G U e \ ([— 1 ,— b] U [b, 1 ]). As before, g e has purely imaginary traces on ((— 1 ,— b) U 
(b, 1 )) n U e that differ by a sign. Moreover, since g e vanishes at e as (z— e) 3 / 2 , g 2 ^ 3 is 
conformal in U e . Furthermore, 


g 2/3 ({x: sgn(e)x < b}n U e ) 

c 

{z : z > 0}, 

,e /3 (((- 1 - b) U (b, 1 )) n U e ) 

c 

{z : z < 0}, 

g 2/3 (A 2± nU e ) 

c 

{z: -sgn(e)arg(z) = ± 27 t/ 3 }. 


where the first property follows from (6.6) and (6.9), the second is a consequence of the 
fact that g e has purely imaginary traces there, and the third is a requirement we impose 
on the arcs A 3 ±. Choosing the branch of g'/ b which is positive on {x : sgn(e)x < b}(T U e , 
one has on ((— 1 ,— b) U (b, 1 )) fl U e that 

(8.17) g 1 + =-sgn(e)ig^ 6 . 
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As in the previous cases, one can verify that the solution of RHP-P e is given by 
(8.18) Pe := E e T 2 ®e (( 3 u/ 2 ) 2 / 3 g 2/3 ) W e , 

where <D_ b := ® and ® b := 03®cr3, W e is defined by (8.12), and 

(8, 9 ) U := MT 2 (p ^ 2 (3n 9e /2)"»/« -L (_^ (e)i 

whose holomorphy can be checked as in the previous cases using (8.17). 


8.4.4 RHP-P±i C 

To map ® into U e , e £ {±icj, set sgn(±ic) = ± and define 

1 rz 


g e (z) := -3 


(hi — h 2 ) (t)dt = -4 log (® (1) (z)/® (2) (z)) / 


z £ U e \ F. Since V is the branch cut for hi and h 2/ the traces of g e on P Hi U e differ 
by a sign. Moreover, since g e vanishes at e as (z— e) 3 / 2 , is conformal in U e - The 
following are conditions we impose on the arcs F, P 2 , and Ti ±: 

ge /j (r 2 nu e ) c {z:z>o}, 

g 2/3 (rnu e ) c {z : z < 0 }, 

ge /j (ri± nu e ) C {z : sgn(e)arg(z) = ± 2 tt/ 3 }. 

Choosing the branch of g'/ 1 ’ which is positive on P 2 P U e , we see that (8.17) holds on 
P P U e with —sgn(e) replaced by sgnl'e). Then the solution of RHP-P e is given by 

(8.20) Pe := E e Ti ®e (( 3 u/ 2 ) 2 / 3 g 2/3 ) W e , 

where <D ic := ® and ®_i c := ff3®(J3, and 


(8.21) 


We := 


:= diag ((® (0) ) n , (® (1) ® (2) ) n/2 v/pl/P 2 , (® (1 '0 (2) ) n/2 v/p 2 /Pl) , 


:= MTi ((p 2 / Pl ( 3 ng e / 2 ) CT 3/6 ^ 


-sgn(e)i 


-sgn(e)i 1 


8.5 Final R-H Problem 

Consider the following Riemann-Hilbert problem (RHP-R): 

(a) R is analytic in C \ Ir, where Lr is the contour shown in Figure 14, and Rf 00) 

(b) R has continuous traces on each side of that satisfy 

I l/VlU‘]l7llVLU'-' _ 

R += R -1 . _i 


= 1 ; 


on n I z , 


(MD rL )j z (MD r 

P e (MD n ) _1 on 3 U e , e £ {± 1 , ±q, ±b, ±ic}, 

where J z was defined in (8.8), while M and D were introduced after (8.9). 



Figure 14. The contour I R : solid lines. The dashed lines represent the relevant 
borders of the domains see Figure 3(a). 
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Let us prove that the jumps of R are uniformly close to I for n e N*. In Section 8.3 we 
have shown that det(N) =, while det(D) = 1 by (2.10) and det(C) = 1 by (2.13). Hence, 
det(M) = 1 and therefore M" 1 = adj(M), where adj(M) is the adjoint matrix of M. Thus, 

(8.22) M(l + 0 (-))M~ 1 =I + MO(-)adj(M) =I + O (0 

uniformly away from {± 1 , ±a, ±b, ±ic} by (7.30). Therefore, it follows from RHP-P e (d) 
that 

(8.23) P e (MD n ) -1 = I + MO( 1 /n)M -1 =I + 0 ( 1 /n) 
uniformly on each 3 U e . Furthermore, we can write 

Jz = I T Pj,ktj,k, j f b, j, k £ { 0 , 1 ,2}, 

where E ,y is the matrix with all zero entries except for the (j + l,k+ 1 )-st one, which is 
1, and pj k is always a combination of pi and P2 (particular values of j, k and the value of 
the entry depend on the arc in question). Thus, 

/®m\ n 

D-JzD- = I+(^yJ Pi,icEyk = I + 0(C^) 

uniformly on Zz (T Zr for some constant Cr > 1 by Theorem 3 (it is a simple examination 
of the five relevant cases). Therefore, we get from (8.22) that 

(8.24) (MD n )j z (MD n ) _1 = I + 0 (C~ n ) 

uniformly on Zz Cl Zr. The relations (8.23) and (8.24) together with [12, Corollary 7.108] 
imply that RHP-R is solvable for all n £ N* large enough and satisfies 

(8.25) R = 1 + 0 ( 1 /n), R(oo) = I, 
uniformly in C, that is, including the boundary values on Zr. 


8.6 Asymptotics of Hermite-Pade Approximants 


Inverting (8.5) and (8.3), we get from (5.4) that 

1 


(8.26) 

where pj = pi; 

(8.27) 


Qn(z) = [Z] n (z) + 


P*(z) 


p(l), x , N Pi (Z) 

R k z = [Z] 12 (z)- 

n Pi z 


^[Z]li4-l (z), Z <E O^-j-, 

0, otherwise, 

—[Z] 13 (z), zeO 0 , 

[Z]13(z), zeOi, 

0, otherwise. 


for z ^ Ft; and 
(8.28) 


*#?(*) = 


0, z e O i 
[Z] 13 (z), otherwise 


+ 


P2lz) 

Pi(z) 


[Z] 12(z), Ot U 0 2 
0, otherwise 


for z f F 2 . 

Let R be the solution of RHP-R. Then 

' MD n in C\U e U e , 


(8.29) 


Z = C" n R 


P e in U e , e G {± 1 , ±a, ±b, ±ic}. 


solves RHP-Z for all n G N, large enough. Denote the first row of R by (1 +u n/ o i> rl j u tV/2 ). 
It follows from (8.25) that 


(8.30) 


|u n ,k(oo)| = 0 and |'U n ,k(z)| = 0 (l/n) uniformly for z £ C, 
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meaning that u n/k (z) is replaced by u^ k (z) when z G Ir. If z ^ U e , e G {± 1 , ±a, ±b, ±ic}, 
then 

2 

(8.31) [Z] lk+ 1 (z) = C n < k ’(z) + C n ^v n j[z)¥V{z), k G { 0 , 1 , 2 }, 

i=o 

according to (4.2), where the functions are replaced by their traces when necessary. 

Formula (8.31) is valid for z G Ui c U U_i c when k = 0 and for z G U e , e G (±b, ± 1 }, 
when k = 1 . Indeed, in these cases P e = MT k+ i ( )W e , see (8.20) and (8.21), (8.11) and 
(8.14), (8.18) and (8.19). As the (k + l)-st column of W e is the same as the (k + l)-st 
column of D n , the (k + 1 )-st column of P e is the same as the (k + 1 )-st column of MD n 
from which the claim follows. 

To estimate the sum on the right-hand side of (8.31), recall that 

(8.32) ^=^,003 0^0-0, 

see (7.27), where 0j0 n/ j0~Q is a rational function on 91 with the divisor 'I) n ,j + 00 !0 - 1 — 
T>n,0 ~ 00 f i 1 . Then (8.30), the maximum modulus principle, and the same normal family 
argument as in (7.31) and (7.32) imply that 

(8.33) |u n/j (z)0 j (z t1c) )0 n/ j(z [1c) )0“ 1 o (z (k) )| = O(l/n;B T1/0 n 9 : t (k) ) 

uniformly in C, where the function on the left-hand side of (8.33) needs to be replaced by 
its traces when z G Ir U 71 ( 391 ^). By combining (8.32) with (8.33) we get that 

(8.34) ^n,j(z)^(z)=¥^(z)O(l/n;D n/0 n 9 t( k )), 3,k G { 0 , 1 , 2 }, 
where O(-) is uniform in C. 

The first relation in (4.4) follows immediately from the first line of (8.26), (8.31), and 

(8.34) . Moreover, (8.31) implies that the second line of (8.26) can be rewritten as 

2 

Qft(z) = C n «J + (Z) +<J_(z)) + C n 2 > n#j (z) «j + (z) +<j_(z)) 

1=0 


forzG Oi + UOi_U O2+ U O2- and z ^ U e , e G {±b, ±q, ± 1 }, where 

z G Oi±, 






T<-W/pt(z), ZGO 


t G {1/2}. 


r=F' 


Clearly, each extends to a holomorphic function in 0 ; + U O U A? by (3.6) (recall 

that V n 1 and ^ n ,2 also satisfy these relations). The first part of (4.5) now follows from 

(8.34). 

Furthermore, (8.27), (8.31), and (8.34) imply the second line of (4.4) outside of Oo U Oi 
for i = 1 . In the spirit of (4.3), define 


$nJ+W = 


/(2) 


ZG Oo, 


Pl/P2)(z)/ ZGOt, 


which is holomorphic in (Oo U Oi U Aq) \ A 2 i by (3.6), and 




/( 2 ), 


zGOi, 


V n/j (z)(Pl/P2)(z), ZGOo, 

which is holomorphic in (Oo U Oi U Aq) \ A-| , again, by (3.6). Then (8.27) and (8.31) imply 
that 

2 


R 


t 3 (z) = c n (w^ 0+ (Z) + <yW_ (Z)) + C n X U n,i (z) (^j + (z) + (z)) 

i=0 
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for z £ (Oo U Oi U Aq) \Fi and z ^ Ll e , e £ {— 1 ,— b, — a,±ic}. Clearly, formula (4.5) is an 
immediate consequence of (8.34)- To finish the proof of (4.4) for i = 1 , let us show that 

^J+ = ^_0(l/n;D n , 0 n>) 

locally uniformly in (i nt (P-| U A01-) UA01-) \ (fl U U ic U Ui c ). Indeed, we have that 


W 

V 


(1) 

n,j+ _ Pi w 2 

(1 ) p 2 W! 

TL.O — 


'( 5 ( 2 )' 

om 


I S P S q n 0 n,j Q i 
1 c(l)c(l)c,(1) 

U n,0 


see (7.27). The claim now follows from Theorem 3 and the normal family argument along 
the lines of (7.3 1 ) and (7.32). This proves (4.4) outside of int|T| U Aqi + ). Analogously we 
can argue that 

$i!J- = ^j + 0(Vn; D n , 0 n IK' 2 )) 

locally uniformly in (int(F-| U Aqi + ) U Aqi + ) \ (Fi U U_i c U Ui c ), which finishes the proof 
of (4.4) for rL 1 1 . The proof of (4.4) and (4.5) for R | 2 * can be completed analogously starting 
with (8.28). 


9 RIEMANN-HILBERT ANALYSIS: CASE II 


9.1 Global Lenses 


Let G1 (u) and G 2 (v) be defined by (8.1). Further, let 0 1 :={z: Re(z) < 0 } be the left half¬ 
plane. We orient the boundary of Oi, say A01 (the imaginary axis), so that Oi lies to the 
left of Aq 1 when the latter is traversed in the positive direction. Denote by Aq 2 a simple 
Jordan curve lying within the right component of H021 containing all the singularities 
°f P2/PI/ see Figure 15 (this is always possible because of Condition 8). We orient A02 



Figure 15. The lens without the circle {|z| = R}, the domains O] (shaded 
regions on the left) and 0 2 (unshaded region), local lenses Ai ± and A 2 ±. 


counter-clockwise and set 0 2 to be the intersection of the exterior domain of A02 and the 
right half-plane (Re(z) > 0 }. Put 


( 9 -i) 


S := Gi (—p/ 2 )G 2 (l/p)Y ■ 


G2(—P2/PI ) G1 (P1/P2) 
f^2( P2/Pi ) 

I 


in Oi, 
in 0 2 , 

in C \ (Oi U 0 2 ) , 


where p := pi (oo)/p2(oo). Put '■= [— 1 , 1 ] U Aqi U Aq 2- Then S solves RHP-S: 

(a) S is analytic in C \I S and lim z _ >oo ±Re(z)>0 S(z)diag (z _2n ,z n ,z n ) = Gi (=Fp/ 2 ); 

(b) S has continuous traces on := Is \ { ± 1 , ±1 /%/!, 0 } that satisfy 


(c) 


J(P 1 , 0 ) 

on 

Ai \{0}, 

J( 0 , p* 2 ) 

on 

A2, 

Gi (P1/P2) 

on 

A01, 

G2(P2/Pi) 

on 

Aq2/ 


where p 2 is defined by (8.4) and A L in (2.8); 

S satisfies RFiP-Y(c) (see Section 5) with [— 1 , 1 ] replaced by Is- 
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If RHP-S is solvable, then so is RHP-Y, and the solutions are connected via (9.1). 

9.2 Local Lenses 

As in Case I, we introduce additional arcs A-\ ± and systems of two arcs A2± as in Figure 15, 
all oriented from left to right. We further denote by the open sets bounded by A[ and 
the arcs Ai±, i £ { 1 , 2 }. Set 

( 9 - 2 ) Z:=SLf in O i± , 

where the matrices are defined by (8.6). Put L z := Ls U Ai + U Ai _ U A2+ U A2—. Then 
it can readily be checked that Z solves RHP-Z: 

(a) Z is analytic in C \ I z and lim z ^ OO/±Re(z)>0 Z(z)diag (z _ 2 n ,z n ,z n ) = G-] (=Fp/2); 

(b) Z has continuous traces on each side of := L z \ { ± 1 , ±1 / \fl, 0 } that satisfy 

' Ji on A? \{ 0 }, i 6 ( 1 , 2 }, 

Jz on ULl( A 0iUA i+ UAi_), 
where J| are defined by (8.7) and 

{ G1 (p i / P2) on A01, 

G2(p2/Pl) on A 02 / 

L t on A i± , i £ { 1 , 2 }; 

(c) Z satisfies RHP-S(c) with replaced by 

As before, if RHP-Z is solvable, then so is RHP-S, and the solutions are connected via 
(9.2). 


Z+ =Z_ 


9.3 Global Parametrix 

Let N = C _n MD n be given by (8.9). Then it is a solution of the following Riemann- 
Hilbert problem (RHP-N): 

(a) N is analytic in C \ [— 1 , 1 ] and lim z _>oo N(z)diag (z _2n , z n ,z n ) = I; 

(b) N has continuous traces on each side of A?, I £ { 1 , 2 }, that satisfy N + — N _ Ji- 

We cannot argue that det(M) = 1 as in Section 8.3 since every entry of M behaves 
like (z=F l/\/ 2 ) -1 / 3 as z -> ±l/-\/ 2 - However, we can construct M 1 explicitly. Denote 
by N = C _n MD n the matrix that satisfies RHP-N as above with pi and pi, replaced by 
1 /pl and 1 /p2■ It follows from (7.29) that the construction (8.9) of the matrices N and N 
is simultaneously applicable or not applicable for each index n. Observe that 


MM 


-1 


MM 


= M 


-1 


(D_/D 


,2n 


M, = M 


-1 


M 


on 7t(cxi) and 7r(( 3 ^) by (2.11) and (2.15). Moreover, 


|mm t J |mm t | = Ji M I = ( 


--1 


JiM, = ( M + J| M_ =1 


--I 


on A? since = J| .It also follows from (3.5) that the entries of M and M have at most 
1/4 and 1 / 3 -root singularities at ±1 and ±1 /\fl, respectively (the functions T tl j and T rli 2 
possess exactly the same behavior around those points as Hence, the product 

MM is holomorphic in the entire complex plane. Since M(oo) = M(oo) = I, we deduce 
that 


I = M M. 


(9-4) 
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In particular, it follows from (3.5) that 
M(z) = (z + 1 

( 9 - 5 ) 


2) 1/3 M± + 0(1), 


M T (z) = (z+l/V 2 ) V3 M± + 0(1), 


as z 


± 1 /V 2 , 


for some constant matrices M± and M ± . Then it follows immediately from (9.4) that 


(9.6) 


M ± M ± = 0 . 


9.4 Local Parametrices 

Again, we need to solve RHP-Z locally, this time around ± 1 , ± 1 /\/ 2 , and 00. The local 
problems RHl’-Pj. 1 are exactly the same as in Case I and therefore their solutions are 
given by (8.i2)-(8.i4). 


9.4.1 RHP-Poo 

Let R > 0 be large enough so that A02 C {|z| < R}. In this section, we are looking for a 
solution to the following Riemann-Hilbert problem (RHP-Poo): 

(a) Poo is holomorphic in {|z| > R) \ Aqi and 

lim C _n Poo(z)diag fz~ 2 n ,z n ,z n ) =G 1 (+p/ 2 ); 

z—>oo,±Re(z)>0 V / 


(b) Poo has continuous traces on each side of A01 D {|z| > R} that satisfy P00+ = 
Poo-Gi (P1/P2); 

(c) Poo = M (i + 0 (n _1 D n uniformly on {|z| = R}. 

Let us show that RHP-Poo, is solved by 

(9-7) Poo := MGt ^(Pi/p 2 )(C 2 /C 1 ) rL u(y^n/2C)'j D n , 


where the constants Cio are defined in (2.13), the function u is given by 

1 2c i f erfc( — V 2 C), Re(C) < 0 , 

2 [ —erfc(\/ 2 C), Re(£) > 0 , 


and C is defined by 


C{z) '■= i^/log (O^ 1 H z )C 2 / 0 ( 2 ) (z)Ci), |z| > R. 

Indeed, £(z) is a conformal function in {|z| > R} that vanishes at infinity by (6.14) (make R 
larger if necessary). Here, we choose the branch of the square root so that z£(z) tends to 
(12) 1 / 4 when z —t 00. Hence, we can deform Aqi in {|z| > R} so that C(Aqi ) C ilR. Thus, 
the right-hand side of (9.7) is holomorphic in {|z| > R) \ Aqi . As it follows from [16, Eq. 
(7.2.2)] that 

lim u(C)=+1/2, 
z—»0,±Re(t)>0 

RHP-Poo(a) is indeed satisfied. To verify RHP-Poo (b), notice that G1 (f + ) = G1 (f_) G1 (f+ — 
f_) and that u + (x) — u_ (x) = e 2x for x C iR. Since 


exp | 2 (n/ 2 K 2 (z)j = (o (2) (z)Ci/O n > (z)C 2 ) 
RHP-Poo, (b) follows. Finally, [16, Eq. (7.12.1)] implies that 


u(C)~ 



L(-n k+1 

k=0 


r(k+ 1/2) (2k+1j 

2 k + 1 r(i/2) 


uniformly in the left and right half-planes. Hence, 


Poo = M (I + (pi/p 2 )(C 2 /C 1 ) n O(l/v / H)E 3/2 ) L) n 
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uniformly on \z\ = R. As |C1 1 = IC2I by Theorem 3, we see that RHP-P 00(c) holds as well. 


9 . 4.2 RHP P.^1 /y/2. 

Denote by U ± 1 := {z : | z 4= 1 / \/ 2 | < n~ 1 / 2 }. In this section, we construct a solution 

to the following Riemann-Hilbert problem: 


(a) P ±l/v / 2 is ana lybc in U ±1 /y/1 \ I z ; 

(b) P ±1 / s/i has continuous traces on each side of L z (T U ± 1 ^ v / 2 that satisfy RHP-Z(b); 

(c) P ±1 /^/2 has the behavior near ±1 /\/ 2 within LI ± ^ ^/ 2 described by RHP-Z(c); 

(d) P ±1 ^y 2 = (^1 + 0(n -1 / 6 )^ MD n uniformly on 9U ±1// ^/2 \ ^-Z- 

In [15] (an alternative approach to asymptotics of multiple orthogonal polynomials 
around cubic branch point was developed in [36]), a 3 x 3 matrix function was constructed 
out of solutions to — Ty'(z) + y (zj =0 that solves RHP-T: 


(a) T is holomorphic in C \ (L+ U I_U (—00,00)), where L± := {£ : Re(£) = ±lm(£)} 

and the positive direction on all the lines goes from the left half-plane to the right 
half-plane; 

(b) T has continuous traces on (L + UL_U (—00,00)) \ { 0 } that satisfy 


t+ = r_ 


and 


T+ = T_ 


T 2 

t 2 

t 3 

t 3 


1 0 

1 1 

0 1 

-1 0 


1 0 

1 1 

0 1 

-1 0 


on L±n{C: Re (0 > 0 }, 
on (0,00), 

on L±n{C: Re (0 < 0 }, 
on (—oo,0); 


(c) Y(C) = O(loglCI) as C-^0; 

(d) Y has the following behavior near 00: 


Y(C;t) = A(C;t) (h-T, (t)C 1/3 +O (V 2/3 )) exp |-|c 2/3 B 2 -t^/ 3 B J 

uniformly for C G C \ (L+ U I_U (—00,00)) and t on bounded sets with 

7-2 


Yi (t) = —- 


T / T 


+ 1 B 


9 


C, 


where 


A(C;x) := \r^ 1/G T 2 (t" 3 / 3 ) < 


B := 


' /_ e 47ii/ 3 1 e 2ra/ 3 N 

1 -1 -1 

y_ e 27ii/ 3 1 e 4ni/3 

e 27ii/ 3 -| e 4ra/ 3 N 

-1 -1 -1 

e 4ra/ 3 ^ e 2ra/ 3 

diag (e 4 ™/ 3 ,1, e 2 ™/ 3 ) , lm(£) > 0, 

diag(e 2 ™ /3 ,l,e 47Ti / 3 ), lm(£) < 0, 


lm(C) > 0, 


lm(C) < 0, 
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and 


C := 


0 

g 4 ra /3 _ g27ii/3 

1 _ e 2 m /3 

e 4 m /3 _ 1 

0 

1 _ g 2 ra /3 

1 _ g47ti/3 

g 4 ra /3 _ g27ii/3 

0 

0 

g 4 ra /3 _ g27ii/3 

e 4 ra /3 _ 1 

1 - e 2 ™ /3 

0 

1 - e 47Ti/3 

e 2 ra /3 _ 1 

g 2 ra /3 _ g47ii/3 

0 


lm(C) < 0. 


Below, we explain how to solve RHP-P 1 using Y, a solution of RHP-P _ 1 -can be 
constructed analogously. 

To carry RHP-T into a neighborhood of 1 / \/ 2 , define 


C(z) := ( \ ^ 


1 /V 2 


1 \“ 1/3 HV 3 ( X ) \ 3/2 
dx 


V2 






where the function H was introduced in (6.16) and we choose the principle branch of 
the square root. Then C is conformal in some neighborhood of 1 / a/2, £(1 /a/2) = 0, and 
C(x) > 0 , x > 1 /V 2 . Further, define 

1 \ 1/3 H _1 / 3 (> 


■=-r~ V3 f 


t(z) := -c 


1/V2 


V2 


x/1 — X 2 


dx. 


It readily follows that x(z) is also conformal in some neighborhood of 1 / a/2, t( 1 /x/2) = 0, 
and that 


(9.8) 


exp <| — -C 2 / 3 (z)B 2 -x(z)C 1/3 (z)B \ = D n (z) 


A+ = A_ 


T 2 

t 3 


on (0,00), 
on (—oo,0). 


by (6.16) and (6.18), and since ® has value 1 at the point of SK whose natural projection is 
l/x/ 2 - Set 

E*(z) := M(z)diag(l, pj -1 (z), pj 1 (z))A _1 (n 3/2 /(z);n 1/2 x(z)). 

It can easily be verified that 

0 r 

,-l Oy 

0 r 

-1 o y 

As the entries of M as well as the entries of 1 can have at most cubic root singularity 
at 1 / a/2, the matrix E* is holomorphic i n 11 ], yy ■ Then 

P*(z) := E*(z)T(n 3/2 C(z);n 1 /2 T(z))diag(l, pi (z), p 2 (z)) 
satisfies RHP-P 1 , (a-c) (we always can adjust so that / maps them into {Re(z) = 

±lm(z)}). It also follows from RHP-T(d) and (9.8) that 

p*(z) = M(Z) (l+Fn(z) +0(u- 2 / 3 )) D u (z) 
as n 4 00 uniformly on cl 11 1 , ^ (recall that t(z) is conformal in LI 1 ^ and vanishes at 

1 /a/2, which implies that n 1 / 2 x(z) remains bounded as n 4 00 and therefore RHP-T(d) 
is applicable), where 

F n (z) := n“ 1/2 C _1/3 (z)diag(l, pf 1 (z),p2 ' 1 (z))Y, (n 1/2 T(z))diag(l, p, (z),p 2 (z)). 
Since T-| (n 1 / 2 x(z)) ~ n 1 / 2 x(z) as z —> 1 /a /2 orn4 00, we can write 


(MfnM-'Kz) = 


r(z)r 1/3 (z) /4 


(z-1/a/2) 2 / 3 


-M+F* M 


iI + o(|z-1/a/2| 1/3 ^ 
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by (9.4) and (9.5), where F* is a matrix of constants (more precisely, it is one matrix of 
constants in {lm(z) > 0} and another one in {lm(z) < 0}, see the definition of Yi) with 
entries whose moduli are uniformly bounded as n —> 00. As the fraction in the above 
equality is holomorphic around 1 /\fl and has value — 3/4 at 1 /\/ 2 , we get that 

(MF n M“ 1 )(z) =-M + F;M++o(|z-l/V 2 | 1/3 ) . 

Since |z — 1 /V 2 | = n -1 / 2 on we have that 

P*(z) = ^1 — M+Fj/M+ + 0 (n -1 / 6 )^ M(z)D n (z) 
uniformly on 911^^. Thus, according to (9.6), RHP-P is solved by 

p 1a /i(z) := (i + m + f;m^p*(z). 


9.5 Final R-H Problem 


The final Riemann-Hilbert problem is RHP-R from Section 8.5 with )z defined in (9.3). 
The same analysis shows that the jump matrices in RHP-R are of order 0 (n~' /6 ) for 



Figure 16. The contour Ir without the circle {|z| = R} (solid lines) and the 
borders of the domains (dashed lines), see Figure 3(b). 


N* 3 n 4 00 and therefore RHP-R is solvable for all n £ N* large enough and satisfies 
(8.25) with n _1 replaced by n _1 / 6 . 


9.6 Asymptotics of Hermite-Pade Approximants 

It can readily be verified that formulae (8.26), (8.27), and (8.28) remain valid. Let R be 
the solution of RHP-R. Then the solution of RHP-Z is given by (8.29), where e £ { ± 
1, ±l/-\/2,00}, for all n £ IN, large enough. The same argument as before shows that 
(8.34) holds in this case as well. 

Since the first column of Poo is the same as the first column of MD n , (8.31) remains 
valid when k = 0 . Thus, the proof of the first parts of (4.4) and (4.5) is exactly the same. 
The second parts of (4.4) and (4.5) are claimed to hold uniformly on each compact subset 
of the respective domains (because contains the point at infinity and Ng l is bounded). 
Hence, given a compact subset, we always can enlarge R in RHP-Poo so that this set is 
contained in {|z| < R}. This way (8.31) is valid on this compact and the remaining part of 
the proof is the same as in Case I. 


10 RIEMANN-HILBERT ANALYSIS: CASE III 

10.1 Global Lenses 

Let Gi (u) and G2 (v) be as in (8.1). Fix a domain, say 0 1, that contains [— 1 , —b] and whose 
boundary Aqi := 30 1 is smooth, lies entirely in Ooi 2 (except for the point where it crosses 
(—b, b)), see Figure 3(c), while crossing the real line at the origin, see Figure 17. In Case 
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Ilia, denote by A02 a smooth Jordan curve lying within the right component of do 21 and 
containing within its interior all the singularities of P2/P1 (recall that P2/P1 = const in 
Case Illb by Condition 8), see Figure 17. We also denote by O2 the intersection of the 
exterior domains of Aqi and Aq 2- We orient Aqi and A02 counter-clockwise and set 



Figure 17. The lens Iz, the domains Oi (shaded region on the left) and Oi (un¬ 
shaded unbounded region), local lenses Ay. and A2±. The curve A02 is present 
only in Case Ilia. 


(10.1) 

Put I s 


S:= G 2 ((P2/Pi)(oo))Y. 


C>2(—P2/P1 )Gi (P1/P2) in Oi, 
G2(-P2/Pl) in 0 2 , 

I in 


C\(Oi U 0 2 ). 


:= [— 1 , 1 ] U Aqi U Aqj. Then, according to (8.2), S solves RPiP-S: 


(a) S is analytic in C \ Is and lim^oo S(z)diag (z 2 n ,z n ,z TL ) = I; 

(b) S has continuous traces on 1 5 := Is \{± 1 , ±a, 0 } that satisfy 


J(P 1 , 0 ) 

on 

\{ 0 }, 

J( 0 , p* 2 ) 

on 


Gi (P1/P2) 

on 

A01, 

G2(P2/Pi) 

on 

Aqi, 


where p{ is defined by (8.4) and A L in (2.8); 

(c) S satisfies RHP-Y(c) (see Section 5) with [— 1 , 1 ] replaced by !$■ 


If RHP-S is solvable, then so is RHP-Y, and the solutions are connected via (10.1). 


10.2 Local Lenses 

As usual, we introduce additional arcs Ai-j- and systems of two arcs A2± as in Figure 17, 
all oriented from left to right. We further denote by Oy the open sets bounded by A[ and 
the arcs Ay-, i £ { 1 , 2 }. Set 

(10.2) Z^SLT 1 in Oi.±, 

where the matrices L, are defined by (8.6). Put Iz := Is U A-| + U A-| _ U A2+ U A2_. Then 
it can readily be checked that Z solves RHP-Z: 

(a) Z is analytic in C \ Iz and lim z _^ O0 Z(z)diag (z 2n ,z _n ,z _rL ) = I; 

(b) Z has continuous traces on each side of I^ := Iz \ {±1 , ±a, ±b, 0 } that satisfy 

Ji on A° \ {0}, i e {1,2}, 

Jz on I^\(A 1 UA 2 ), 


Z+ =z_ 
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(10.3) 


where are defined by (8.7) and 
r 


Jz := 


J(P1,0) 

on 

(-1/ 

Gi (P 1 /P 2 ) 

on 

A 01 , 

£>2(P2/p1 ) 

on 

Ao2/ 

Li 

on 

^i±/ 


(c) Z satisfies RHP-S(c) with Ig replaced by Iz- 


As before, if RHP-Z is solvable, then so is RHP-S, and the solutions are connected via 

(10.2). 


10.3 Global Parametrix 

Let N = C~ n MD n be given by (8.9). Then det(N) = 1 and it is a solution of the following 
Riemann-Hilbert problem (RHP-N): 

(a) N is analytic in C \ (Ai U A2) and lim z _ s . 0O N(z)diag (z” 2rL ,z n , z n ) = I; 

(b) N has continuous traces on each side of A° that satisfy N + = N_ J L , i £ { 1 , 2 }. 

10.4 Local Parametrices 

Again, as in the previous cases, we need to solve RHP-P e for e £ {± 1 ,±a, ±b}. In fact, 
these local problems are exactly the same as in Case I. Thus, their solutions were con¬ 
structed in Section 8.4. 

10.5 Final R-H Problem 

Once more, the final Riemann-Hilbert problem is RHP-R from Section 8.5 with Jz defined 
in (10.3) and e £ {±q, ±b, ± 1 }. Exactly the same analysis shows that the jump matrices in 



Figure 18. The contour I R (solid lines) and the relevant borders of the domains 
CLjk (dashed lines), see Figure 3(c,d). 


RHP-R are of order 0 ( 1 /n) for n £ IN* and therefore RHP-R is solvable for all n £ IN* 
large enough and satisfies (8.25). 

10.6 Asymptotics of Hermite-Pade Approximants 

As in Case I, one can verify that the formulae (8.26), (8.27), (8.28), (8.31), and (8.34) remain 
valid in this case as well. So the proof of (4.4) and (4.5) proceeds exactly along the same 
lines as in Case I. 






io.6 ASYMPTOTICS OF HERMITE-PADE APPROXIMANTS 


REFERENCES 


[1] A. Angelesco. Sur deux extensions des fractions continues algebriques. Comptes Rendus de 1 'Academie des 
Sciences, Paris, 18:262-265, 1919. 3, 24 

[2] A.I. Aptekarev. Asymptotics of simultaneously orthogonal polynomials in the Angelesco case. Mat. Sb., 
i36(i78)(i):56-84, 1988. English transl. in Math. USSR Sb. 64, 1989. 3, 4, 24 

[3] A.I. Aptekarev. Strong asymptotics of multiple orthogonal polynomials for Nikishin systems. Mat. Sb., 

1 9°(5) : 3 _ 44/ x 999- English transl. in Sb. Math. 1999. 3, 4 

[4] A.I. Aptekarev. Asymptotics of Hermite-Pade approximants for a pair of functions with branch points. Dokl. 
Akad. Nauk, 422(4)1443-445, 2008. English trans. in Dokl. Math. 78(21:717-719, 2008. 10 

[5] A.I. Aptekarev, P.M. Bleher, and A.B.J. Kuijlaars. Large n limit of Gaussian random matrices with external 
source. Part II. Comm. Math. Phys., 259:367-389, 2005. 4 

[6] A.I. Aptekarev, V.A. Kalyagin, V.G. Lysov, and D.N. Toulyakov. Equilibrium of vector potentials and uni- 
formization of the algebraic curves of genus 0. /. Comput. Appl. Math., 233(602-616), 2009. 24 

[7] A.I. Aptekarev, A.B.J. Kuijlaars, and W. Van Assche. Asymptotics of Hermite-Pade rational approximants 
for two analytic functions with separated pairs of branch points (case of genus o). Int. Math. Res. Papers, 
2008, rpmoo7, 128 pp. 4 

[8] A.I. Aptekarev and V.G. Lysov. Systems of Markov functions generated by graphs and the asymptotics of 
their Hermite-Pad approximants. Mat. Sb., 20i(2):29~78, 2010. English transl. in Sb. Math., 201:183-234, 2010. 

3 

[9] A.I. Aptekarev and H. Stahl. Asymptotics of Hermite-Pade polynomials. In A.A. Gonchar and E.B. Saff, 
editors. Progress in Approximation Theory, pages 127-167, Berlin, 1992. Springer-Verlag. 4 

[10] A.I. Aptekarev, D.N. Toulyakov, and W. Van Assche. Hyperelliptic uniformization of algebraic curves of the 
third order. J. Comput. Appl. Math., 284:38-49, 2015. 25 

[11] A.I. Aptekarev, D.N. Toulyakov, and M.L. Yattselev. On uniformization of a certain cubic algebraic curve of 
genus 2. manuscript. 20 

[12] P. Deift. Orthogonal Polynomials and Random Matrices: a Riemann-Hilbert Approach, volume 3 of Courant Lec¬ 
tures in Mathematics. Amer. Math. Soc., Providence, RI, 2000. 41 

[13] P. Deift, T. Kriecherbauer, K.T.-R. McLaughlin, S. Venakides, and X. Zhou. Strong asymptotics for polynomi¬ 
als orthogonal with respect to varying exponential weights. Comm. Pure Appl. Math., 52(i2):i49i-i552, 1999. 

39 

[14] P. Deift and X. Zhou. A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics for 
the mKdV equation. Ann. of Math., 137:295-370, 1993. 14 

[15] K. Deschout and A.B.J. Kuijlaars. Double scaling limit for modified Jacobi-Angelesco polynomials. In 
P. Branden, M. Passare, and M. Putinar, editors. Notions of Positivity and the Geometry of Polynomials, Trends 
in Mathematics, pages 115-161, Basel, 2011. Springer. 46 

[16] F.W.J. Olver et al. editors. NIST digital library of mathematical functions, http://dlmf.nist.gov. 38, 45 

[17] U. Fidalgo Prieto and G. Lopez Lagomasino. Nikishin systems are perfect. Constr. Approx., j^(f}:2()y-^6, 
2011. 3 

[18] A.S. Fokas, A.R. Its, and A.V. Kitaev. Discrete Panleve equations and their appearance in quantum gravity. 
Comm. Math. Phys., 142(2):313~344, 1991. 13 

[19] A.S. Fokas, A.R. Its, and A.V. Kitaev. The isomonodromy approach to matrix models in 2D quantum gravi¬ 
tation. Comm. Math. Phys., 147(2):395~430, 1992. 13 

[20] F.D. Gakhov. Boundary Value Problems. Dover Publications, Inc., New York, 1990. 27 

[21] A.A. Gonchar and E.A. Rakhmanov. On convergence of simultaneous Pade approximants for systems of 
functions of Markov type. Trudy Mat. Inst. Steklov, 157:31-48, 1981. English transl. in Proc. Steklov Inst. Math. 
15Z 1983- 3/ 24 

[22] A.A. Gonchar, E.A. Rakhmanov, and V.N. Sorokin. Hermite-Pad6 approximants for systems of Markov-type 
functions. Mat. Sb., i 88(5):33-58, 1997. English transl. in Math. USSR Sbornik i88(5):67i-696, 1997. 3, 24 

[23] C. Hermite. Sur la fonction exponentielle. C. R. Acad. Sci. Paris, 77:18-24, 74-79, 226-233, 285-293, 1873. 2 

[24] V.A. Kalyagin. On a class of polynomials defined by two orthogonality relations. Mat. Sb., no(4):6o9-627, 
1979. 3, 24 

[25] A.B.J. Kuijlaars, K.T.-R. McLaughlin, W. Van Assche, and M. Vanlessen. The Riemann-Hilbert approach to 
strong asymptotics for orthogonal polynomials on [—1,1]. Adv. Math., i88(2):337-398, 2004. 37 

[26] A.A. Markov. Deux demonstrations de la convergence de certaines fractions continues. Acta Math., 19:93- 
104, 1895. 3 

[27] E.M. Nikishin. A system of Markov functions. Vestnik Moskovskogo Universiteta Seriya 1, Matematika 
Meklianika, 34(4):6o-63, 1979. Translated in Moscow University Mathematics Bulletin 34(4), 63-66, 1979. 3, 24 

[28] E.M. Nikishin. Simultaneous Pade approximants. Mat. Sb., H3(i55)(4):499-5i9, 1980. 3 

[29] J. Nuttall. Asymptotics of diagonal Hermite-Pade polynomials. J. Approx. Theory, /[2(f):2<)()-^86, 1984. 4 

[30] H. Pade. Sur la representation approchee d'une fonction par des fractions rationnelles. Ann. Sci Ecole Norm. 
Sup., 9(3):3-93, 1892. 2 

[31] I.I. Privalov. Boundary Properties of Analytic Functions. GITTL, Moscow, 1950. German transl., VEB Deutscher 
Verlag Wiss., Berlin, 1956. 22 


I 10.6 ASYMPTOTICS OF HERMITE-PADE APPROXIMANTS 


[32] E.A. Rakhmanov. On the asymptotics of Hermite-Pade polynomials for two Markov-type functions. Mat. 
Sb., 202(1)1133-140, 2011. English transl. In Sb. Math., 202(i):i2y-ij 4, 2011. 3 

[33] H. Stahl. Asymptotics of Hermite-Pade polynomials and related approximants. A summary of results. Man¬ 
uscript, 79 pages. 4 

[34] H. Stahl. Asymptotics of Hermite-Pade polynomials and related convergence results - a summary of results. 

In Nonlinear numerical methods and rational approximation, volume 43 of Math. Appl., pages 23-53. Reidel, 
Dordrecht, 1987. 4 

[35] K. Strebel. Quadratic Differentials, volume 5 of Ergebnisse der Mathematik und Hirer Grenzgebiete (3). Springer- 
Verlag, Berlin, 1984. 21 

[36] D.N. Tulyakov. Difference schemes with power-growth bases perturbed by the spectral parameter. Mat. Sb., 
20o(5):i29-i58, 2009. English transl. in Sb. Math., 200(51:753-781, 2009. 46 

[37] W. Van Assche, J.S. Geronimo, and A.B.J. Kuijlaars. Riemann-Hilbert problems for multiple orthogonal 
polynomials. In Special functions 2000: current perspective and future directions, number 30 in NATO Sci. Ser. II 
Math. Phys. Chem., pages 23-59, Dordrecht, 2001. Kluwer Acad. Publ. 13 

[38] E.I. Zverovich. Boundary value problems in the theory of analytic functions in Holder classes on Riemann 
surfaces. Russian Math. Surveys, 26(1)1117-192, 1971. 27 

Keldysh Institute of Applied Mathematics, Russian Academy of Science, Moscow, Russian Federa¬ 
tion 

E-mail address: aptekaa@keldysh.ru 

Department of Mathematics, KU Leuven, Celestijnenlaan 200B box 2400, BE-3001 Leuven, Belgium 

E-mail address: walter@wis.kuleuven.be 

Department of Mathematical Sciences, Indiana University-Purdue University Indianapolis, 402 North 

Blackford Street, Indianapolis, IN 46202, USA 

E-mail address: maxyatts@math.iupui.edu 


